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ABSTRACT 

We probe the geometry of magnetospheric accretion in classical T Tauri stars by modeling red 
absorption at He I A10830 via scattering of the stellar and veiling continua. Under the assumptions 
that the accretion flow is an azimuthally symmetric dipolc and helium is sufficiently optically thick that 
all incident l-/^m radiation is scattered, we illustrate the sensitivity of He I A10830 red absorption to 
both the size of the magnetosphere and the filling factor of the hot accretion shock. We compare model 
profiles to those observed in 21 CTTS with subcontinuum redshifted absorption at He I A10830 and 
find that about half of the stars have red absorptions and l-/jm veilings that are consistent with dipole 
flows of moderate width with accretion shock filling factors matching the size of the magnetospheric 
footpoints. However, the remaining 50% of the profiles, with a combination of broad, deep absorption 
and low 1-^m veiling, require very wide flows where magnetic footpoints are distributed over 10-20% 
of the stellar surface but accretion shock filling factors are < 1%. We model these profiles by invoking 
large magnetospheres dilutely filled with accreting gas, leaving the disk over a range of radii in many 
narrow "streamlets" that fill only a small fraction of the entire infall region. In some cases accreting 
streamlets need to originate in the disk between several i?* and at least the corotation radius. A few 
stars have such deep absorption at velocities > 0.5 V^ sc that flows near the star with less curvature 
than a dipole trajectory seem to be required. 

Subject headings: accretion, accretion disks — planetary systems: protoplanetary disks — scattering 
— stars: formation — stars: pre-main-sequence 



1. INTRODUCTION 

The classical T Tauri stars (CTTS) are optically 
revealed, low-mass, pre-main-sequence stars that ac- 
crete material from a circumstellar disk and have a 
well-defined connectio n between accretion and outflow 
(jHartigan et al.l [19951 HEG hereafter). Accretion from 
the disk to the star is thought to be guided by 
the stellar magnetosphere, where a sufficiently strong 
magnetic field truncates the disk at several stellar 
radii and material follows field lines t hat direct it to 
the stellar surface at high latitudes dGhosh fc Lama 
1975 iKonigll [19911 ICoifier Cameron fc Campbelll 11993 
Shu et aLlll994T h Magnetospheric accretion controls the 
star-disk interaction, and an improved understanding of 
this process will shed light on outstanding issues in the 
innermost 10 R* of CTTS systems, such as the regula- 
tion of stellar angular momentum and the launching of 
the inner wind. 

Leading diagnostics of CTTS accretion include the op- 
tical/UV continuu m excess and the pr ofiles of permit- 
ted emission lines (jBouvier et al.ll2007a| ). Kinematic ev- 
idence for infalling gas in CTTS began with the dis- 
covery that some CTTS show inverse P Cygni struc- 
ture in upper Balmer lin es extending to velocities of 
several hundred km s _1 (jWalkerl Il972f ). Later, more 
sensitive surveys found that redshifted absorption com- 
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ponents are relatively common in some l ines, especially 
in the upper Balmer an d Paschen series ( Edward s et al 
119941: lAlencar fc Basril 120001 : iFolha fc Emersonl .2001 
Although redshifted absorption extending to several hun- 
dred km s" 1 clarifies that material accretes in free-fall 
from at least several stellar radii, it has been the success 
of radiative transfer modeling of line formation in mag- 
netospheric ac cretion flows in a key s eries of papers cul- 
minating with iMuzerolle et al.l (|2001h that has provided 
the strongest underpinning for this phenomenon. Under 
the assumption of an aligned, axisymmetric dipole, the 
models have had reasonably good success in reproducing 
the general morphology of hydrogen profiles and emission 
fluxes in some stars. The complementary assessment of 
accretion rates follows from interpreting the SED of the 
optical/UV excess, which has been successfully modeled 
for wavelengths shortward of 0.5 /im as arising in a hot 
accretion shock, where accreting material impacts the 
stellar surfa ce after free-fall along f unnel flows coupled 
to th e disk (|Johns-Krull et alJl2000HCalvet fc Gullbrind 
119981 CG hereafter). To match the observations, the 
accretion shock filling factor is less than 1% in most 
cases but can climb to 10% in a few of the most ac- 
tive accretors. The derived accretion rates range from 
10~ 10 to 10~ 6 Mp. yy- 1 with a median of 10~ 8 M Q yr" 1 
(jCalvet et al]l2000h . 

Additionally, Zeeman broadening of unpolarized CTTS 
photospheric lines indic ates mean su r face fi eld strengths 
in the range 1-3 kG ()Johns-Krulll 120071 ). sufficiently 
strong to induce disk truncation and drive funnel flows. 
However, these strong surface fields are not predomi- 
nantly dipolar, as photospheric lines show only weak net 
circular polarization implyi ng dipole components an or- 
der of magnitude smaller (jValenti fc Johns-K rull 2004; 



I Yang et aI1l2007t ). Nevertheless, an extended dipole com- 
ponent is inferred for the accretion flow, since the same 
authors find significant circular polarization in the nar- 
row component of the He I A5876 emission line, thought 
to be forme d in the accretion sh ock at the base of the 
funnel flow (jBeristain et al.ll200ih . 

The evidence for magnetospheric accretion is thus com- 
pelling; however, the topology of the magnetosphere, the 
geometry of the accretion flow, and the disk truncation 
radius remain topics of considerable investigation, since 
their configuration impacts processes for angular momen- 
tum regulation and wind launching. One form of an- 
gular momentum regulation, known as disk locking, in- 
vokes a spin-up torque from accreting material just inside 
the corotation radius balan ced by a spin-down to rque at 
larger radii (|Collier Cameron fe Campbell I1993T). This 
approa ch has been questioned bv lMatt fc Pudritzl (|2005l 
2OO8a0), who instead suggest accretion-powered stellar 
winds as a more likely means for stellar spindown, which 
must occur simultaneously with ma gnetospheric accre - 
tion from the disk (see also ISautv fc Tsinganosj :1 994). 
Alternatively the X-wind model (|Shu et al.l [1994) orig- 
inally featured a narrow annulus of star-disk coupling 
close to the corotation radius, where closed field lines 
develop funnel flows and open field lines drive a centrifu- 
gal outflow that carries away angular momentum from 
accreting material, thus inhibiting stellar spin-up. The 
flexibility of this model to maintain its basic properties 
in the face of complex magnetospheric accretion geome- 
tries h as recently been demonstrated by lMohantv fc Shul 
( 2008) . If the stellar and disk fields are parallel, then an 
intermittent outflow can develop via a Reconnection X- 
wind, which removes an gular momentum fro m the star 
as well as the inner disk (jFerreira et al.ll2006t ). 

Evidence for non-aligned fields coupled with complex 
accretion geometries is mounting, coming from a vari- 
ety of recent studies. Using time-resolved spectropo- 
larimetry of the m ildly accreting CTTS V2129 Oph, 
iDonati et al.l (|2007f ) used Zeeman detections from both 
photosphcric features and emission lines from the accre- 
tion shock to construct a Doppler tomographic map of 
the magnetic topology on the stellar surface. The dom- 
inant field on the star is a misaligned octupole, and ac- 
cretion is confined largely to a high-latitude spot, cov- 
ering < 5% of the stellar surface. These authors also 
attempt to reconstruct the 3D field geometry out to the 
disk interaction region and suggest that the large-scale 
field funneling accreting material from the disk is more 
complex than a simple dipole. However, such extrap- 
olation techniques, while tantalizing, must be applied 
with caution at present, since there are numerous uncer- 
taint ies in the reconstruction process (|Mohantv fc Shul 
2008). Numerical magnetohydrodynamic (MHD) simu- 
lations of star-disk interactions demonstrate that non- 
axisymmetric funnel flows arise if the stellar field is 
tipped by only a few degrees relative to the rotation 
axis, breaking into t wo discrete streams un der stable ac- 
cretion conditions (Ro manova et al.l l2003f l. If the ac- 
cretion rate is sufficiently high, accretion is predicted 
to proceed throug h equatorial "tongues " that can push 
apart field lines ([Romanova et al.ll2008h . Similarly, ac- 
cretion spots can appear at a wide range of latitudes, 
including equatorial belts. Such behavior has been ob- 
served in models featuring accretion along quadrupo- 



lar a s well as dipolar field lines (|Long et alj 120071 
l2008h . accretion alo ng field lines extrapo lated from sur- 
face magnetograms ([Gregory et al.l loops'), and accretion 
mediated by a dynamo-generated disk magnetic field 
(|von Rekowski fc Brandenburg||2006f h 

Observational signatures for misaligned dipoles are be- 
ing explored in radiative transfer models for hydrogen 
line formation in funnel flows, with subsequent predic- 
tions for rotationally modula ted profile variations. Ini- 
tially (SynfingtOTLetlal] (2005) presented radiative trans- 
fer models of hydrogen lines featuring curtains of accre- 
tion covering a limited extent in azimuth in geometries 
consistent with aligned dipoles. Their model profiles ex- 
hibit certain characteristics of the observed line profile 
variability, such as rotationally modulated line strengths 
and the appearance of red absorption components at 
certain phases and inclinations, but the predicted level 
of variability is highe r than observed. More recently, 
iKurosawa et al.1 (|2008h applied a radiative transfer code 
for H line for mation to the 3-D output fro m the MHD 
simulations of Ro manova et al.l (|2003l . l2004f ). which pre- 
scribe the geometry, density, and velocity of two-armed 
accretion streams that result from dipoles misaligned 
with the rotation axis by angles ranging from 10° to 
90°. Applying temperatures similar to those from the 
iMuzerolle et al.l (|2001h axisymmetric models, they were 
able to reproduce some of the trends in continuum and 
profile variability from models. One of the larger dis- 
crepancies in comparing the model profiles to observed 
ones is that the model profiles for Paschen and Brackett 
lines are a factor of two narrower t han the mean value ob- 
served by [lHha~&jimersoil (|200l|) . The problem is likely 
more complex than simply finding another line broaden- 
ing mechanism, since Paschen 7 has recently been shown 
to have line widths that are correlated with the l-//m 
continuum excess, in the sense that the narrowest lines 
are fo und among objects w ith the lowest disk accretion 
rates (lEdwards et al.|[20fM EFHK hereafter). Evidence 
that some of the hydrogen emission is formed in the ac- 
cretion shock rather than the funnel flow is now clearly 
demonstrated from the discov ery of circular polar ization 
in the core of B aimer lines (|Donati et al.l 120081 ) . The 
implication is that hydrogren lines are not necessarily 
a definitive means for probing the properties of funnel 
flows, so additional probes are desirable. 

In this paper, we explore a different means of diag- 
nosing the geometry of the accretion flow, making use 
of the redshifted subcontinuum absorption in the 2p 3 P° 
— > 2s 3 S transition of neutral helium (A10830), recently 
demonstrated to be a very sensitive probe of both out- 
flowing gas in the inner wind and infalling gas in the fun- 
nel flow due to its frequent display of blue and red ab- 
sorptions (EFHK). A subcontinuum absorption feature 
is a more tell-tale diagnostic of a kinematic flow than an 
emission profile, since its position (blue or red) indicates 
the direction of the flow, its width indicates the range of 
line-of-sight velocities in the flow, and its depth at a par- 
ticular velocity indicates the fraction of the continuum 
(stellar plus veiling) occulted by material moving at that 
velocity. In the case of a red absorption, the absorption 
depth signals the fraction of the stellar surface covered 
by the funnel flow at each velocity, making it an effective 
probe of the CTTS accretion geometry. 

This is the third in a series of papers about l-/mi di- 



agnostics of accretion and outflow in CTTS. The first, 
EFHK, presented l-/^m spectra from 38 CTTS including 
He I A10830 profiles, P7 profiles, and measurements of 
the continuu m excess "veil i ng" in the 1-/Ltm region. The 
second paper, Kwan et al. (2007, KEF hereafter), mod- 
eled blueshifted absorption components at He I A 10830, 
which appear in about three quarters of the sample, and 
found that while some stars have winds best explained 
as arising from the inner disk, others require an out- 
flow moving radially away from the star in an accretion- 
powered "stellar" wind. In this paper, we analyze red- 
shifted subcontinuum absorption at He I A10830 in 21 
CTTS that present red absorption in at least one obser- 
vation. The following section describes the sample se- 
lection, data acquisition, and data reduction. Section 3 
presents the data and discusses variability. In Section 4, 
we present model scattering profiles that arise in a dipo- 
lar flow geometry, show that they explain only a fraction 
of the red absorptions, and explore modifications to a 
dipolar flow that better explain the remaining observa- 
tions. Discussion and conclusions follow in Sections 5 
and 6. 

2. SAMPLE AND DATA REDUCTION 

In this paper we focus on the 21 of 38 CTTS included in 
EFHK that display redshifted subcontinuum absorption 
at He I A10830 at least once in a multi-epoch observing 
program with Keck NIRSPEC. It includes spectra pre- 
sented in EFHK acquired in November 2001 and Novem- 
ber 2002, when 8 of the 38 CTTS were observed twice 
and 1 on three occasions. It also includes 33 additional 
spectra of 24 objects from that study, taken in 2005, 
2006, and 2007. In EFHK, 19 out of 38 CTTS showed 
red absorption at He I A10830. In the subsequent ob- 
serving runs, He I A10830 red absorption was seen in 2 
additional stars. Thus, the 21 of 38 CTTS (55%) that 
have shown subcontinuum red helium absorption in at 
least one spectrum of 81 acquired between 2001 and 2007 
form the sample for this paper. Among the 21 stars with 
helium red absorption, 12 were observed more than once, 
with 6 observed twice, 1 observed three times, 4 observed 
four times, and 1 observed six times. The EFHK sam- 
ple was assembled to span the full range of mass accre- 
tion rates observed for CTTS, from less than 10~ 9 to 
~ 10~ 6 M Q yr" 1 , with a median rate of 10~ 8 M Q yr _1 . 
Most of them are from the Taurus- Auriga star-forming 
region and have spectral types of K7 to M0. The subset 
of 21 stars that are the focus of this paper are identi- 
fied in Table [TJ along with their spectral types, masses, 
radii, rotation periods, median veilings ry at 0.57 /zm, 
and mass accretion rates from the literature. We have in- 
cluded ry only for the 29 sources in common with HEG, 
obtained more than a decade earlier. We also list the 
number of observations for each star. 

As in the 2006 paper, the additional sp ectra were ac- 
quired with NIRSPEC on Keck II f;McLean et all 1 19981) 
using the Nl filter (Y band), which covers the range 
0.95 to 1.12 ^m at a resolution R = 25,000 (AT/ = 
12 km s~ l ) . The echelle order of primary interest extends 
from 1.081 to 1.096 ^m and contains both He I A10830 
and P7. Spectra from the 2005-06 season were acquired 
by G. Blake (13 December 2005) and D. Stark (13 Jan- 
uary 2006). Those from November and December 2006 
were obtained by L. Hillenbrand, W. Fischer, S. Edwards, 



and C. Sharon. Finally, Hillenbrand obtained two addi- 
tional spectra of TW Hya in December 2007. Data reduc- 
tion, including wavelength calibration and spatial recti- 
fication, extraction of one-dimensional spectra from the 
images, and removal of telluric emission and absorption 
features, is discussed in EFHK. While we used an IRAF 
script to reduce the EFHK data, we used the IDL pack- 
age REDSPEC by S. S. Kim, L. Prato, and I. McLean to 
reduce data acquired in Fall 2006 and later. EFHK also 
describes the procedure for measuring photospheric lines 
to determine the l-/um veiling ry, defined as the ratio 
of excess flu x to photospheric flux near the He I A10830 
line (see also lHartigan et aLlll989l ). After the veilings are 
determined, a non-accreting template that has been ar- 
tifically veiled to match the observed CTTS is subtracted 
from each target spectrum. This removes photospheric 
absorption lines from the He I A10830 and P7 regions, 
which allows for a more accurate definition of the remain- 
ing structure in each of these two lines. 

We augmented the spectral templates from those of 
EFHK, resulting in a reassessment of the l-/im veil- 
ing for one object. Recent determinations of the spec- 
tral type of BM And in th e V ba nd range from G 8 
(jGuenther fc Hessmanl [1991 to K5 (jMora et al.ll2001h . 
In EFHK we used an early K star to deveil the 2002 
spectrum of BM And. However, using our new grid of 
templates acquired in Fall 2006, we found that the G8 
dwarf HD 75935 provides a better match to the photo- 
sphere of BM And. Deveiling the 2002 spectrum of BM 
And with this template yields a veiling of 0.4, in contrast 
to the value of 0.1 reported in EFHK. We adopt the more 
appropriate veilings for BM And in this work, ry = 0.4 
in 2002 and ry = 0.5 in 2006. The veiling determina- 
tions for all the other objects from EFHK are unaffected 
by our extended grid of templates. 

We use the stellar mass, radius, and rotation period 
to calculate the escape velocity and the star-disk corota- 
tion radius for each star, which are included in Table [1] 
and will be used in later analysis of the accretion ge- 
ometry. We carefully surveyed the literature to acquire 
the most up-to-date es timates of spectral types and stel- 
lar luminosities , using iKenvon fe Hartmannl (|1995f ) and 
iGullbring et al.l (|1998l ) in most cases. The luminosity 
of YY Ori (HEG) was updated to reflect the latest es - 
timate of the distance to Orion (|Menten et al.l 12007) . 
which is 10% less than the earlier value. Spectral types 
were converted into effective temp eratures using the scale 
from iHillenbrand fc White! (J2004I ). and stellar radii fol- 
low directly from application of the Stefan-Boltzmann 
law to the effective temperatures and luminositie s. Stel - 
lar masses are then derived from the iSiess et all (2000) 
pre-main-sequence tracks, available online. Since the es- 
cape velocity will be an important parameter in com- 
paring observed to model profiles, we have given some 
thought to its accuracy. With a dependency on M/R, 
the largest source of uncertainty in calculating the es- 
cape velocity is the uncertainty in T c g, since temperature 
strongly affects both the mass and radius determination, 
while luminosity only weakly influences the radius es- 
timate. We assess that the typical error in the escape 
velocity is ~ 20%. For the 12 stars with rotation peri- 
ods in the literature, we calculate corotation radii with 
a typical error of 20%, provided the photometric period 
is equivalent to the rotation period. 



TABLE 1 

21 OF 38 CTTS WITH SUBCONTINUUM RED ABSORPTION AT He I A10830 



Object 


Sp Type 


M, 


R* 


Vcsc 


Prot 


Rco 


(rv) log M aC c Ref 


N obs 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


AA Tau . . 


K7 


0.70 


1.75 


390 


8.22 


8.7 


0.32 


-8.5 


12,9,5,9 


4 


BM And . . 


G8 


2.03 


3.02 


510 








>-9 


15,15,7 


2 


CI Tau . . . 


K7 


0.70 


1.94 


370 






0.47 


-6.8 


12,12,10 


1 


CY Tau . . 


Ml 


0.43 


1.70 


310 


7.5 


7.2 


1.20 


-8.1 


12,9,4,9 


3 


DK Tau . . 


K7 


0.69 


2.51 


320 


8.4 


6.1 


0.49 


-7.4 


12,9,3,9 


4 


DN Tau . . 


MO 


0.52 


2.15 


300 


6.0 


5.2 


0.08 


-8.5 


12,9,2,9 


2 


DR Tau . . 


K7 


0.69 


2.75 


310 


9.0 


5.9 


9.60 


-5.1 


12,10,4,10 


4 


DS Tau . . . 


K5 


1.09 


1.30 


570 






0.96 


-7.9 


12,9,9 


1 


FP Tau... 


Ml 


0.21 


2.00 


200 






0.15 


-7.7 


12,12,10 


1 


GI Tau . . . 


K6 


0.93 


1.74 


450 


7.2 


8.8 


0.24 


-8.0 


12,9,17,9 


2 


GK Tau . . 


K7 


0.69 


2.16 


350 


4.65 


1.8 


0.23 


-8.2 


12,9,3,9 


2 


HK Tau . . 


M0.5 


0.45 


1.65 


320 






1.10 


-6.5 


12,12,10 


1 


LkCa 8 . . . 


MO 


0.53 


1.48 


370 


3.25 


5.1 


0.15 


-9.1 


12,9,3,9 


2 


RW Aur B 


K5 


0.96 


1.09 


580 








-8.8 


19,19,19 


1 


SU Aur . . . 


G2 


2.02 


3.27 


490 


1.7 


2.3 




-8.0 


12,12,6,8 


1 


TW Hya.. 


K7 


0.75 


1.04 


520 


2.80 


7.3 




-9.3 


18,18,13,14 


6 


UY Aur . . 


MO 


0.54 


1.30 


100 






0.40 


-7.6 


11,11,11 


1 


UZ Tau E 


Ml 


0.43 


1.39 


310 






0.73 


-8.7 


19,19,19 


1 


UZ Tau W 


M2 


0.33 


1.88 


260 








-8.0 


11,11,11 


1 


V836 Tau. 


K7 


0.71 


1.43 


440 


7.0 


9.6 


0.05 


-8.2 


12,12,16,10 


2 


YY Ori . . . 


K7 


0.68 


3.00 


290 


7.58 


4.8 


1.80 


-5.5 


10,10,1,10 


1 
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Note. — Col. 2: Spectral type; Col. 3: Stellar mass in Mq; Col. 4: Stellar radius in Rq; Col. 5: Stellar escape 
velocity in km s _1 , calculated from columns 3, 4; Col. 6: Rotation period in days; Col. 7; Corotation radius in R t , 
calculated from columns 3, 4, 6; Col. 8: Median veiling at 5700 A from HEG; Col. 9: Logarithm of the mass accretion 
rate in Mq yr — 1 ; Col. 10: References for the spectral type, stellar luminosity (to determine M* and R*), rotation 
rate (where available), and mass accretion rate; Col. 11: Number of spectra acquired with NIRSPEC. 



Three of the 21 objects are known members of binary 
pairs resolved in our spectra where we have observed only 
the primary: DK Tau A, HK Tau A, and UY Aur A. For 
another system, RW Aur, we have resolved spectra of 
RW Aur A and RW Aur B, but only the latter shows red 
absorption at He I A10830 and thus qualifies as part of 
the sample for this study. There is conflicting evidence in 
the literature on whether RW Aur B has a close compan- 
ion at an angu lar separation of 0.12" with a K-band flux 
ratio of 0.024 (|Ghez et al.llJ99l iCorreia et al1l200rl . In 
our spectra we see the lines of only one object consistent 
with a K5 spectral type, and we call this RW Aur B. 
An additional two objects are unresolved binaries: UZ 
Tau E and UZ Tau W. For these, we also see lines from 
only one star and attribute the l-/um continuum and line 
profiles to the primary. 

3. EMPIRICAL RESULTS 

In this study we concentrate on the red absorption seen 
at He I A10830 as a probe of the accreting gas, ignor- 
ing the blue absorptions that arise from disk and stellar 
winds (KEF). For each spectrum of the 21 stars that 
show subcontinuum redshifted He I A10830 absorption 
at least once, Table [5] lists the HJD of observation, the 
l-/im veiling, and measurements of the red absorption. 
In this section we first compare the veilings of the stars 
in this study to those from the ensemble of 38 CTTS in 
EFHK, and then we present the profiles and kinematic 
data for the reference sample, consisting of the single ob- 
servation of each star with the deepest red absorption 
(identified with an asterisk in Table [5]). Next we demon- 
strate that the propensity for He I A10830 to absorb all 



impinging l-/im photons provides a means of estimating 
the origination radius in the disk for infalling gas and 
the filling factor of accreting material immediately be- 
fore the accretion shock. We conclude the section with a 
discussion of profile and veiling variability. 

3.1. Veiling and Redshifted Absorption 

Our additional observations beyond those in EFHK 
confirm the result reported therein, that subcontinuum 
redshifted absorption is more prevalent in CTTS with 
low veiling. We illustrate this in Figure [TJ where the 
equivalent width of the redshifted He I A10830 absorption 
below the continuum is plotted against both the simulta- 
neous l-/xm veiling ry and the non-simultaneous optical 
veiling r v for the 38 CTTS in EFHK. The 21 CTTS that 
are the focus of this study, showing redshifted absorption 
at least once among 46 spectra, are each identified by 
name. The remaining 17 CTTS that have not yet been 
seen to show redshifted absorption among 35 spectra ac- 
quired to date appear as symbols (but can be identified 
from EFHK). All points in the figure are averages for 
objects with multiple observations taken between 2001- 
2007. 

We have included the optical veiling measurements in 
Figure Q] because it is the excess emission at optical and 
shorter wavelengths that is associated with luminosity 
from accretion shocks and is the basis for deriving disk 
accretion rates. Note that the range of veilings is dif- 
ferent at each of the two wavelengths, with maxima of 
ry = 9.6 and ry = 2, and that although all CTTS show 
detectable veiling in the optical (HEG), 7/38 have no 



TABLE 2 
Veilings and Measurements of He I A10830 Subcontinuum Red Absorption 









W x 


■'-'max 


V C 


FWQM 


^blue 


Vi-cd 


Object 


HJD 


ry 


(A) 


(%) 


(km s" 1 ) 


(km s — *) 


(km s — 1 ) 


(km s — ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


AA Tau . . 


605.0 


0.2 


0.5 


14 


185 


120 


110 


250 




606.9 


0.1 


1.9 


12 


110 


170 


10 


250 




1718.0* 


0.0 


1.5 


61 


90 


310 


-10 


310 




2069.0 


0.1 


0.9 


21 


50 


70 


a 


80 


BM And . . 


604.8 


0.4 


2,3 


28 


115 


260 


-20 


290 




2068.7* 


0.5 


2.9 


10 


90 


270 


-10 


280 


CI Tau . . . 


605.9 


0.2 


1.3 


17 


110 


290 


10 


310 


CY Tau . . 


606.8 


0.1 


1.0 


27 


110 


120 


80 


230 




1718.0* 


0.0 


1.1 


37 


110 


100 


80 


210 




2068.8 


0.2 


0.0 












DK Tau . . 


604.9 


0.5 


2.1 


28 


115 


290 


20 


330 




606.9 


0.5 


1.9 


37 


80 


210 


a 


280 




1748.9* 


0.0 


3.1 


10 


150 


320 


-20 


310 




2068.9 


0,1 


1.9 


34 


95 


290 


a 


310 


DN Tau . . 


606.0 


0.0 


1.3 


30 


115 


170 


60 


250 




1718.0* 


0.0 


1.2 


33 


150 


110 


70 


260 


DR Tau . . 


605.0 
606.0 


2.0 
2.0 


0.0 
0.0 














606.9* 


2.0 


0.7 


11 


235 


160 


150 


320 




2069.1 


3.5 


0.0 












DS Tau... 


605.9 


0,1 


1.1 


18 


205 


210 


90 


310 


FP Tau... 


605.0 


0.1 


0.5 


17 


50 


120 


a 


120 


GI Tau . . . 


606.0 


0.1 


3.1 


17 


160 


230 


50 


3.30 




2069.8* 


0.0 


3.3 


52 


180 


210 


50 


350 


GK Tau . . 


606.0 


0.3 


0.0 














2069.8* 


0.1 


0.5 


11 


160 


110 


50 


220 


HK Tau . . 


606.1 


0.4 


0.5 


17 


75 


100 


30 


110 


LkCa 8 . . . 


604.9* 


0.05 


1.1 


32 


160 


160 


70 


280 




2068.9 


0.1 


1.2 


21 


125 


190 


10 


250 


RW Aur B 


605.1 


0.1 


2.8 


13 


160 


230 


50 


3.30 


SU Aur . . . 


607.0 


0.0 


1.6 


35 


50 


180 


-50 


150 


TW Hya.. 


605.2 
606.1 


0.0 
0.0 


0.0 
0.0 














1718.1* 


0.1 


1.1 


32 


245 


170 


170 


370 




2069.1 


0.1 


0.9 


17 


255 


170 


170 


350 




2452.1 


0.0 


0.8 


17 


2.30 


190 


150 


350 




2453.1 


0.1 


0.7 


11 


210 


190 


160 


3.30 


UY Aur . . 


605.0 


0,1 


0,1 


11 


160 


110 


110 


220 




607.0 


0,1 


0,1 


12 


160 


110 


100 


220 




1718.1* 


0.2 


0.7 


20 


160 


130 


90 


210 




2069.9 


0,3 


0.5 


13 


160 


120 


90 


230 


UZ Tau E 


605.9 


0,3 


0.2 


8 


185 


60 


150 


210 


UZ Tau W 


605.9 


0.1 


0.6 


18 


75 


110 


20 


170 


V836 Tau. 


606.0* 
1749.0 


0.0 
0.0 


1.7 
0.0 


35 


160 


170 


80 


300 


YY Ori . . . 


607.1 


0,1 


2.1 


37 


225 


210 


110 


390 



Note. — Col. 2: Heliocentric Julian Date (2,452,000 +); for multiple observations an asterisk indicates membership 
in the reference sample; Col. 3: Veiling at one micron; Col. 4: Equivalent width of red absorption below the continuum; 
Col. 5: Percentage of the continuum absorbed at the deepest point of the profile; Col. 6: Centroid of the red absorption; 
Col. 7: Width at one quarter of red absorption minimum; Col. 8: Minimum velocity of red absorption; Col. 9: 
Maximum velocity of red absorption. 
a The true minimum velocity is obscured by central absorption; we assume Vbi ue = 0. 



detectable veiling at 1 ^m. (All the CTTS with r Y = 
do show He I A10830 and P7 in emission, differentiat- 
ing them from WTTS.) Despite these differences, objects 
with low ry have low ry, and objects with high ry have 
high ry , clarifying that the l-fim veiling is a rough proxy 
for optical/UV veiling and disk accretion rate, and that 
variations in the disk accretion rate are relatively mod- 
est over timescales of a decade. The proportionality be- 
tween ry and disk accretion rate is further corroborated 
by the excellent correlation between ry and the equiv- 
alent width of P7 emission (EFHK), assuming that the 
equivalent widt h of P7, like that of P/3 , is correlated with 
accretion rate (Muzcroll e et al.l ll998: Fo lha fc Emerson! 
I2001L iNatta et al.ll2004h ~ 



The prime message from Figure [T] is that when the 
veiling is high, ry > 0.5 or ry > 2, red absorption at 
He I A10830 is rare. Although the number of observa- 
tions of each of the 38 stars ranges from 1 to 6, we note 
that out of a total of 25 observations of the 9 objects 
with ry > 0.5, only once, in one of four observations of 
the highest-veiling object DR Tau, did a weak redshifted 
absorption appear. In contrast, out of the 56 total spec- 
tra of the 29 objects with ry < 0.5 or ry < 2, redshifted 
absorption is detected in 37 spectra. Even with our non- 
uniform sampling of individual objects, it is clear that 
the frequency of redshifted absorption in CTTS with the 
highest veilings, seen in only 4% (1/25) of the total spec- 
tra of 9 objects, is significantly lower than that in CTTS 
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Fig. 1. — Equivalent width of red absorption at He I A10830 
versus veiling for CTTS from EFHK. The top panel shows the re- 
lation for the simultaneously measured l-/jm veiling ry for all 38 
CTTS, using averages for stars with multiple observations. The 21 
stars featured here show helium red absorption in at least one ob- 
servation and are labeled with abbreviations of their names. The 
17 that have not been observed to show red absorption are iden- 
tified with plus signs. The bottom panel plots the same equiva- 
lent width data versus the average of optical veiling measurements 
(A = 5700 A) from HEG, obtained a decade before the NIRSPEC 
campaign, which exist for 29 of the EFHK stars. In this and future 
scatter plots, points that would otherwise overlap are slightly offset 
for clarity. 

with more modest veilings, where red absorption is seen 
in 66% (37/56) of the total spectra of 29 objects. Spec- 
tral variability for these objects will be discussed in Sec- 
tion 3.4. 

3.2. Line Profiles and Subcontinuum Absorption 

He I A 10830 profiles for the 21 CTTS that have shown 
subcontinuum redshifted absorption at least once are 
presented in Figured] This set of profiles is for the refer- 
ence sample, and the profiles are ordered by their simul- 
taneous 1-fim veiling. The part of the profile we identify 
as the red absorption component is delineated in Fig- 
ure [5] by shading. As noted above, the reference sample 
contains the profile with the deepest red absorption at 
He I A 10830 for each star (in contrast to the reference 
sample from EFHK that emphasized blueshifted absorp- 
tion from winds). The full set of profiles observed for 
the 12 stars with multiple spectra appears in Section [3.4l 
where we discuss variability. The reference sample will 
be used in all subsequent analysis unless we are explicitly 
considering profile or veiling variations. 

The uniqueness of the He I A10830 line in the study 
of accreting gas is immediately seen by comparing it to 
the P7 profile found in the same NIRSPEC order. This 
comparison is made in Figure [3J which zooms in on the 
red side of the profile for each of the 21 objects in the 
reference sample, sorted now by the equivalent width of 
the He I A10830 red absorption. (We have ignored the 
blue half of the line in order to draw attention to the 



red absorption; full P7 profiles can be found in EFHK.) 
Only 5 of 21 (24%) stars show red absorption at both 
He I A10830 and P7, and when seen, it is considerably 
weaker at P7. Specifically, the maximum depth of red 
absorption seen at P7 is 21% of the continuum, compared 
to 61% for A10830, and the maximum equivalent width 
is 1.1 A, versus 4.5 A at A10830. Surprisingly, the 3 
stars with the strongest helium absorption (AA Tau, Gl 
Tau, and DK Tau) show no absorption at P7, while the 
3 stars with the strongest P7 absorption (TW Hya, BM 
And, and YY Ori) have intermediate helium absorptions, 
although their P7 absorptions do share similar velocity 
structures with their helium absorptions. In models for H 
line formation in magnetospheric accretion scenarios, in- 
verse P Cygni absorption is seen only when the accretion 
rate is favorable and the line of sight is directed toward 
the hot continuum in the accretion shock. In contrast, 
our data indicate a much wider range of formation con- 
ditions for red absorption at He I A10830, which offers 
a unique probe of the infalling gas projected in front of 
the stellar surface by absorbing continuum photons from 
both the star and the accretion shock. 

Measured parameters of the He I A10830 red absorp- 
tion, i.e., the section of the profile shaded in Figure [51 
are listed for each observation in Table [2] with an as- 
terisk identifying the spectrum in the reference sample 
for stars with multiple observations. Parameters include 
the equivalent width W\, the depth of maximum pene- 
tration into the continuum D max , the centroid velocity 
Vc, and the width measured at one quarter of the ab- 
sorption minimum, FWQM. We also tabulate velocities 
at the blueward and redward edges of the absorption, 
Vbiuc and Vi-ed- In most cases, Vbiuc is easily identified 
as the location where emission sharply transitions to red 
absorption. On the other hand, the gradual return to 
the continuum at the high-velocity end makes V Te d less 
straightforward to measure. In order to have a uniform 
definition for all stars, we conservatively define V%d as 
the velocity where the absorption reaches 95% of the con- 
tinuum level, with the consequence that it is somewhat 
smaller than the extreme infall velocity. Histograms il- 
lustrating the diversity of these parameters appear in 
Figured The equivalent widths range from 0.2 to 4.5 A, 
maximum penetrations into the continuum range from 
8% to 61%, centroids range from 50 to 255 km s _1 , and 
FWQM range from 60 to 320 km s -1 . In many stars the 
absorptions begin near the stellar rest velocity, so the 
FWQM reflects the true width of the absorbing veloci- 
ties. In others Vbiuc is well redward of the rest velocity 
(e.g., DR Tau and TW Hya) due to the presence of he- 
lium emission that is likely from another region, such as 
the wind, that is filling in the red absorption and reduc- 
ing its magnitude. 



3.3. Maximum Infall Velocities 

Without assuming any particular infall geometry, we 
can estimate the outer extent of an accretion flow by 
comparing the most redward velocity in an absorption 
profile with the stellar escape velocity. A particle under- 
going ballistic infall from a distance R toward a star of 
mass M* and radius R„ has a free-fall speed at a distance 
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Fig. 2. — The reference sample of residual He I A10830 profiles for the 21 CTTS with subcontinuum rcdshiftcd absorption (shaded), 
ordered by decreasing l-/im veiling ry ■ Since the veiled photosphcric contribution has been subtracted, the continuum corresponds to zero 
on the flux axis, and total absorption of the continuum corresponds to -1. Velocities are relative to the stellar photosphere, and the spectra 
are plotted with three-pixel binning. 
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of the red half of He I A10830 (left) and P7 (right) profiles from the reference sample, arranged top to bottom 
A10830 red absorption equivalent width. Subcontinuum absorption is shaded in both lines. 
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Fig. 4. — Histograms of red absorption profile measurements for 
the reference sample. 
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(1) 



where T4 SC is the escape velocity from the surface of the 
star. The largest infall velocity achieved in the funnel 
flow, immediately before impact, is thus set by the max- 
imum distance where infalling gas leaves the disk, i? max , 

1 /2 

i.e., V max = V csc (1 - R*/R max ) ■ We can then use the 
velocity of the red edge of the He I A10830 absorption, 
V rc d, as an indicator for V max and hence determine -R max . 
The depth of the absorption near T4 e d will then indicate 
the filling factor of infalling gas immediately before it im- 
pacts the photosphere. However, because of projection 
effects and the conservative assumption for measuring 
Vrcd, this gives a lower limit to the actual Vmax, and thus 
the corresponding inferred maximum distance for infall, 
Rrcdi will be a lower limit to the actual value of -R max . 
Thus V max > Vrod and with i? max in units of R* , we have 



Rb 



> (1 - V r 2 ed /V e 2 sc ) 



R Te 



(2) 



In Table [3] we list V re d, V, sc , and their ratio, followed 
by the implied R YO d- Figure [5] shows the locations of the 
21 observations of the reference sample in the (V-ed, V, sc ) 
space as well as dotted lines corresponding to V rc d/Ve SC of 
0.94, 0.87, 0.71, and 0.30, or R Ted = 8, 4, 2, and 1.1 R* 
respectively. The average value of i? rod is 2.9 i?», and 
the median is 1.9 i?*, where we adopt R lc d > 8-R* for 
the 3 stars with Vred > Vc SC . Of these 3 outliers, YY 
Ori has V re d/V csc = 1.4, which indicates an error in the 
stellar parameters, while the other two, DR Tau and DK 
Tau, have V-ed slightly larger than Vo SC but within the 
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Fig. 5. — Velocity at the red edge of the He I A10830 absorption 
(^rcd) versus the stellar escape velocity (Vise) for the 21 profiles of 
the reference sample. The solid line represents V rcc i = V esC} while 
the dotted lines represent the ratios of V re( j to Vc SC for ballistic 
infall from 8, 4, 2, and 1.1 ij*, assuming V rc d = Vmax- 

estimated 20% uncertainty. If we used a less conservative 
estimate for V c d as the outermost redward velocity (see 
Section 3.2), at a penetration depth of 2% rather than 
5% of the continuum, then the median _R rc d increases to 
2.9 R*. 

In Figure [S] we compare R Ic d to R co for the 12 stars 
with published rotation periods (listed in Table [lj. For 
this group of stars, the median ratio of R lc d to R co is 0.4, 
which increases to 0.5 for the less conservative estimate 
of the maximum redward velocity and can increase fur- 
ther when projection effects are considered. This is well 
inside the corotation radius for most stars, but as will 
be apparent from model profiles in Section 4, for some 
viewing angles projection effects can result in a signif- 
icant underestimate of -R max as determined from i? re d- 
Three stars show R Icd /R co > 1.3 (YY Ori, DR Tau, and 
DK Tau) , indicating that infall originates close to or pos- 
sibly outside corotation, unless the error in R co is larger 
than the typical 20% uncertainty. 

The projected area of the accretion flow immediately 
above the accretion shock can be equated to the depth 
of the absorption at the highest velocity in the red 
absorption profile. This estimate will be most reliable 
for objects where both (1) V re d is near Vc S c, indicating 
that projection effects are not significantly altering our 
determination of the velocity just before impact, and 
(2) the l-/im veiling is near zero, indicating that the red 
absorption is solely due to scattering of stellar photons 
so the absorption depth at each velocity is the minimum 
percentage of the stellar surface that is obscured at 
that velocity. In the reference sample, three stars with 
Vcd > 0.8 Vcsc and no 1-fim veiling have significant 
absorption depths at 0.9 V rcd : AA Tau (14%), DK Tau 
(10%), and GI Tau (10%). Their implied projected 
areas of material moving close to the escape velocity are 
an order of magnitude larger than the accretion shock 
filling factors estimated from their optical continuum 
excesses (CG). Although accretion flows do "funnel" into 
narrow columns as they arrive at magnetic footpoints 
on the stellar surface, the coverage fraction of the flow 
in typical dipole flow geometries diminishes by less than 
50% from ~ 1.2 i?* to i?* as the speed increases from 
~ 0.9 Vred to Ved, not by an order of magnitude. This 
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TABLE 3 
Maximum Infall Distances and Corotation Radii 



Object 


Vred 


V'csc 


Vrod/Vcsc 


fired a 


i?co 


^Lrcd/^co 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


AA Tau . . 


310 


390 


0.79 


2.7 


8.7 


0.3 


BM And . . 


280 


510 


0.55 


1.1 


[6.3] 


[0.2] 


CI Tau . . . 


310 


370 


0.84 


3,1 


[6.3] 


[0.5] 


CY Tau . . 


240 


310 


0.77 


2.5 


7.2 


0.4 


DK Tau . . 


310 


320 


1.06 


[>8] 


6.1 


[>1.3] 


DN Tau . . 


260 


300 


0.87 


4.0 


5.2 


0.8 


DR Tau . . 


320 


310 


1.03 


[>8] 


5.9 


[>1.4] 


DS Tau... 


310 


570 


0.60 


1.6 


[6.3] 


[0.3] 


FP Tau... 


120 


200 


0.60 


1.6 


[6.3] 


[0.3] 


GI Tau . . . 


350 


150 


0.78 


2.5 


8.8 


0.3 


GK Tau . . 


220 


350 


0.63 


1.7 


1.8 


0.1 


HK Tau . . 


140 


320 


0.44 


1.2 


[6.3] 


[0.2] 


LkCa 8 . . . 


280 


370 


0.76 


2.3 


5.1 


0.5 


RW Aur B 


330 


580 


0.57 


1.5 


[6.3] 


[0.2] 


SU Aur . . . 


150 


190 


0.31 


1.1 


2.3 


0.5 


TW Hya.. 


370 


520 


0.71 


2.0 


7.3 


0.3 


UY Aur . . 


210 


100 


0.60 


1.6 


[6.3] 


[0.3] 


UZ Tau E 


210 


340 


0.62 


1.6 


[6.3] 


[0.3] 


UZ Tau W 


170 


260 


0.65 


1.7 


[6.3] 


[0.3] 


V836 Tau. 


300 


110 


0.68 


1.9 


9.6 


0.2 


YY Ori . . . 


390 


290 


1.34 


[>8] 


4.8 


[>l-7] 



Note. — For the reference sample only. Col. 2: Maximum velocity of He I A10830 red absorption (km s ); Col. 3: 
Escape velocity from the stellar surface (km s — 1 ); Col. 5: Lower limit to maximum distance of infalling material (-R*); 
Col. 6: Corotation radius (-R*). 

a Brackets indicate an assumed R TB d > 8 R* since V YBC i > Vesc- b Brackets indicate that, since P ro t and thus R co are 
unknown, R co is set to the mean of the known values. c Brackets indicate uncertainty due to one of the preceding 
conditions. 

3.4. Variability 
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stars with known P ro t- Sort is by R rcd /R co , which decreases from 
top to bottom. Lower limits to R Ie d are indicated by greater-than 
symbols. 



hints that in at least several stars, a conventional dipolar 
accretion flow will have difficulty reconciling small shock 
filling factors with deep, broad red absorptions. 



The He I A10830 profiles for the 12 CTTS in this study 
with multiple spectra are shown in Figure [7J where for 
each star the full set of observed profiles is superposed 
and the range of simultaneous veilings is indicated. Since 
the time intervals are randomly distributed, ranging from 
days to years, only very general statements about vari- 
ability can be made. Three categories of variability are 
seen: (1) five objects always show redshifted absorption 
with little variation in the absorption morphology (BM 
And, UY Aur, LkCa 8, DN Tau, and GI Tau); (2) two 
objects always show red absorption, but the profile mor- 
phology changes dramatically (AA Tau and DK Tau); 
and (3) five objects have no redshifted absorption at one 
epoch, but do show it at another epoch (DR Tau, GK 
Tau, TW Hya, CY Tau, and V836 Tau). 

The 5 stars observed at least four times (AA Tau, DK 
Tau, DR Tau, TW Hya and UY Aur) can be exam- 
ined to see if there is a relation between veiling and the 
He I A10830 red absorption. The red absorption equiv- 
alent width is plotted against ry for each observation 
of these five stars in Figure El For the star with lit- 
tle change in the morphology of its red absorption (UY 
Aur), the veiling varies by a factor of 2. For the two 
stars where redshifted absorption is always present but 
changes dramatically (AA Tau and DK Tau), the ab- 
sorption is strongest when the veiling is lowest (i.e., not 
detected) . For the two stars where redshifted absorptions 
come and go (DR Tau and TW Hya), there is no relation 
between veiling and the strength of the absorption. 

Each of the 5 stars with at least 4 observations was ob- 
served on at least two nights of a three-night run in 2002, 
providing a look at short-term variability and the possi- 
ble role of rotation. Data points from this run appear 
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Fig. 7.— Residual He I A10830 profiles of the twelve CTTS that were 
absorption in at least one observation. The reference sample spectra are 
appears in each box. 

in Figure [5] as open symbols, and the pair of asterisks 
for TW Hya are from a second set of two consecutive 
nights 5 years hence in 2007. The only objects to show 
much variation over a time scale of days are DR Tau, 
which showed weak red absorption only on the last of 
three consecutive nights in 2002, and AA Tau, which we 
now discuss further. 

The variability of AA Tau at optical wavelengths has 
been thoroughly examined in the context of rotational 
modulation from a mi saligned magnetosphere in t eractin g 
with the inner disk (|Bouvier et alj I1999L [2001 l2007bl ). 
Briefly, the system is close to edge-on, with an inclina- 
tion angle of 75°, and the rotation period is 8.22 days. 
Phase zero corresponds to the epoch of maximum V-band 
flux, while phase 0.5 is characterized by a reduction in 
V-band flux due to occultation of the star by a warped 
disk. Accretion diagnostics are strongest near phase 0.5, 
with redshifted absorption appearing at Ha and H/3 be- 
tween phase 0.39 and phase 0.52 accompanied by a rise in 



observed more than once and that show subcontinuum redshifted 
shown with heavier lines, and the range of observed l-/jm veilings 

the optical veiling (measured between 5400 and 6700 A) 
from 0.2 at phase zero to between 0.4 and 0.7 during the 
occultation phase. 

To see whether our l-/im data of AA Tau are consis- 
tent with this picture, we adopt the 8.22-d ay rotation pe- 
riod, a ssign phase 0.51 to HJD 2,453,308 (jBouvier et al.1 
l2007bf l. and convert our observation dates to rotation 
phases. Figure [5] shows the He I A10830 and P7 profiles 
and veilings for each of our 4 observations, correspond- 
ing to projected phases from ~ to 0.4. The figure 
also plots the equivalent width of He I A10830 red ab- 
sorption against the derived phase, where each point is 
roughly aligned with its corresponding profiles. While 
He I A10830 red absorption appears at all phases, it is 
weakest near phase zero and increases steadily to phase 
0.39. The velocity at the red absorption edge (Ked) also 
varies, increasing from 250 km s _1 near phase zero to 
310 km s _1 at phase 0.39. In contrast, P7 shows red 
absorption only once, close to phase zero, when the 1- 
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SCATTERING MODELS AND COMPARISON TO 
OBSERVATIONS 
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Fig. 8. — Equivalent width of red absorption at He I A10830 
versus the 1-fim veiling for stars with at least four observations 
and at least one helium profile with subcontinuum red absorption. 
Observations on contiguous days are represented by open points 
for 2002 and asterisks for 2007. 
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Fig. 9. — Relation between red absorption and rotational phase 
for 4 spectra of AA Tau (phased from Bouvicr ct al. 20072). The 
He I A10830 and P7 profiles of AA Tau are shown above the cor- 
responding equivalent width of the red absorption at He I A10830 
for each phase. Profiles are labeled with the simultaneous l-/im 
veiling, and their velocity axes run from -500 to 500 km s _1 . For 
the helium line, the flux axis runs from -1 to 1, while for P7, it 
runs from -0.5 to 0.5 to elucidate the morphology of the weaker 
profiles. 



/im veiling is also highest. If the phasing from the Bou- 
vier epoch is accurate, then the deepest and widest red 
absorption at He I A10830 occurs at the phase associ- 
ated with maximum accretion effects in the optical when 
the line of sight pierces the disk warp and the accretion 
shock. However, this would then mean that the l-/im 
veiling and the P7 red absorption are out of phase with 
respect to optical veilings and profiles. Whether or not 
this phase projection is accurate, the profile sequence for 
He I A10830 and P7 provides another illustration of the 
very different kinds of information about the accretion 
flow that can be inferred from these two lines. Clearly, 
time-monitoring studies at 1 /im will be revealing! 



Radiative transfer models of hydrogen 
mg in the accretion flow 



line s aris- 
dMuzerolle et~aT1 120011 



ISvmington et alj|2005t iKurosawa et al.ll2006t |20M ) have 
been successful in reproducing general characteristics of 
these lines in some stars. However, the models assume all 
of the hydrogen emission arises in the funnel flow, they 
depend on an assumed temperature in the flow that is 
not well understood, and they are limited in their ability 
to constrain the accretion geometry. In this section, we 
take a new approach to understanding CTTS accretion 
flows, modeling the scattering of continuum photons by 
He I XI 0830 in the inf ailing gas. The lower level of the 
He I A10830 transition (2s 3 S) is 21 eV above the ground 
state, restricting its formation to regions near the star 
where the ionizing photon flux is high. Further, since 
the only permitted transition downward from the upper 
level (2p 3 P°) is emission of a A10830 photon, we model 
this line as resonance scattering. 

We first lay out the assumptions of our model, in which 
the accretion geometry is the commonly adopted dipolar 
flow, a geometrically flat disk is truncated by the inner- 
most field lines, and all accreting field lines terminate in 
an accretion shock of uniform temperature at the stel- 
lar photosphere that generates a continuum excess ob- 
servable as veiling. We then compare profiles generated 
from these models to the observed He I A10830 red ab- 
sorption profiles. The basic dipolar flow is found lacking 
in a significant number of objects, so we then explore 
modifications to this geometry that better explain these 
observations. 

4.1. Basic Dipolar Flow 

We first consider an axisymmetric dipolar field in 
which the stellar magnetic and rotational axes are aligned 
and an opaque accretion disk extends from an initial ra- 
dius Ri to infinity in the equatorial plane. The outline of 
the overall structure of accretion from the disk to the star 
is completely specified by two parameters, although there 
is some flexibility in which two parameters we choose. 
One pair is Ri and Rf, which indicate the range in ra- 
dial distance from the star, i.e., Ri < R < Rf, over 
which the dipolar field lines that participate in accretion 
are distributed over the disk. Alternatively, we can spec- 
ify 6i and 9f, which mark the range in polar angle, i.e., 
Of <6 < Oi and (% — 6i) < 9 < (it — Of), where the same 
field lines are distributed at the stellar surface. The rela- 
tion between the two pairs is apparent from the dipolar 
field structure, which stipulates that 



Rij = R*/ sin 2 Oi j. 



(3) 



A third pair of parameters is F and Rq, where F = 
cos Of — cos Oi is the fraction of the stellar surface out- 
lined by the above field lines, and R = i?* / sin 2 , with 
cos Of — cos 0q = cos #0 — cos Oi , marks the approximate 
median radius of where the accretion flow originates in 
the disk. The median field line originating at Ro will 
thus correspond to a median polar angle on the star of 
. We find this pair of parameters to be instructive, 
since if the whole geometric structure is fully occupied 
by accreting gas, then F will equal /, the filling factor 
of the accretion shock. In this subsection, we consider 
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Model Magnetosphehic Geometry Parameters 
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Fig. 10. — Schematic representations of accretion geometries 
used in scattering calculations for dipoles with F = f. The star is 
black, the accretion flows are gray, and the disk is the solid line in 
the equatorial plane. Each column shows a different Ro (#o), and 
each row shows a different / with a corresponding ry We note 
that in the extreme case of Ro = 8-R* and / = 0.1, the accreting 
field lines thread the disk out to 35 R, , so the shading in the lower 
right panel extends far beyond the figure boundary. 



F = f and thus use / to indicate both the fraction of 
the stellar surface outlined by the overall flow structure 
and the filling factor of shocked gas at the terminus of 
accreting field lines, as in previous work by others. 

The modeled values of Rq and / are chosen to sam- 
ple the full range of plausible accretion flow sizes and 
filling factors. The values of Rq, taken to be 2, 4, and 
8 R*, are consistent with the understanding that the ac- 
cretion flow arises nea r the star-disk corotation radiu s 
(|Ghosh fe Lambl [19781 iKonigll [l99H iShu et all fl99l . 
The values of /, taken to be 0.01, 0.05, and 0.1, cover the 
range found from shock models of CG. The upper sec- 
tion of Table 0] lists the 9 modeled combinations of Rq 
(0o), / (= F), an d the associated ranges (Ri,Rf) over 
which material leaves the disk. These configurations are 
visualized in Figure [TU] (The lower half of the table lists 
cases with F =/= f, which will be explored beginning in 
Section 4.2.) Three cases correspond closely to geome- 
tries used in previous m odels of magnetospheric accretion 
(jMuzerolle et al.ll200JI) . The case with R = 4.R, and 
/ = 0.01 approximates their SN (small/narrow) case, 
the case with Rq = 4i?* and / = 0.05 approximates 
their SW (small/ wide) case, and the case with R = 8i?* 
and / = 0.01 approximates their LW (large/wide) case. 

For each model the veiling r\ from the associated ac- 
cretion shock, defined as the ratio of the continuum ex- 
cess flux F v to the stellar flux F*, is determined by the 
blackbody temperatures of the star and hot gas from the 
shock, the magnitude of /, and the viewing angle. In 
all cases we assume a T* = 4000 K blackbody for the 
stellar continuum and a T v — 8000 K blackbody for the 
continuum from the accretion shock. This is a typical 
value found by CG from continuum excesses shortward 
of 0.5 /xm, although values as high as 10,000 K or as 
low as 6000 K are sometimes indicated. The veiling at 
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Note. — Col. 1: Fiducial disk coupling radius (i?»); 
Col. 2: Stellar impact angle in degrees from the pole; Col. 3: 
Fraction of the star over which the full range of magneto- 
spheric footpoints is distributed; Col. 4: Filling factor of 
accretion shocks on the stellar surface; Col. 5: Approximate 
1-fim veiling (eq. [2]); Col. 6: Innermost radius at which 
accreting material leaves the disk, also the disk truncation 
radius (-R*); Col. 7: Outermost radius at which accreting 
material leaves the disk (-R*). Ri and Rf follow directly 
from Rq and F. 



wavelength A is 



r A = tt 



IT (T v ) ( f 



If (T*) VI -/ 



(4) 



where the approximate equality arises from setting the 
ratio of the projected areas of the two continua perpen- 
dicular to the line of sight, which depends on viewing 
angle, to simply //(l — /). Over the full range of view- 
ing angle, the observed r\ for the same / can vary by a 
factor of a few (see Section 4.1.1). 

The approximate value for ry from equation ^ , with- 
out including the effect of viewing angle, is identified in 
Figure [10] for each of the 3 values of /. With our as- 
sumed temperatures, the ratio of the blackbody inten- 
sities from the veiling continuum and the photosphere 
Ij* (T v )/I b x b (T*) is 24.5 at A = 5700 A and 6.3 at 
A = 1.08 /im, so that for a typical observed / = 0.01, the 
approximate veilings at these wavelengths are ry = 0.25 
and ry = 0.06. The corresponding ratio of ry fry ~ 4 is 
preserved for all / and is independent of viewing angle. 

The velocity of the flow has contributions from both 
free-fall and rotation. The free-fall speed at a distance r 
from the star along a field line threading the disk at R 
is given by equation ([1]). Since the gas follows the field 
lines, the velocity vector takes the form 



v // 



-Vff 



3(7 1 / 2 (l-( Z ) 1 / 2 / o±(2-3g)z 



(5) 



(4 - 3<z)V2 

(|Calvet fc Hartmann|[l992l : lHartmann et al.lll994D . Here 
q = sin 6, and (p, <p, z) are unit vectors in the cylindrical 
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coordinate system. Above the equatorial plane, the plus 
sign applies, while below the equatorial plane, the minus 
sign applies so that the flow is always from the disk to 
the star. 

For the rotational component of the flow, the magne- 
tosphere is assumed to rotate rigidly with velocity 



v* -5-0, 
it* 



(6) 



where v* is the rotation speed of the star at its equator, 
assumed here to be 0.05 V = c , or 15 km s _1 when VL 



300 km S -1 , a typical value for TTS (jRebull et al.ll2004 
and references therein), and p is the cylindrical radial 
distance of a point from the rotation axis. Since the 
rotational motion is for the most part transverse to the 
line of sight for the absorbing gas seen projected in front 
of the star, it has a very small effect on the absorption 
part of the line profile. 

The flow scatters continuum photons, which arise from 
the star and the accretion-heated photosphere. To max- 
imize the red absorption, the A10830 transition in the 
accreting flow is assumed to be optically thick. A rectan- 
gular line absorption profile with a 10 km s _1 half-width 
is adopted to account for thermal and turbulent broaden- 
ing. Thus, if a particular ray from a point on the stellar 
surface intersects the accreting flow such that the projec- 
tion of the gas velocity along the ray extends from v m { n to 
ti max , then continuum photons from (w m in — 10 km s -1 ) 
to (wmax + 10 km s _1 ) are scattered. Because sponta- 
neous emission is the dominant de-excitation route of the 
A10830 upper state (2p 3 P°) in comparison with other 
decay, collisional, or ionization processes, the photon ab- 
sorption and subsequent re-emission is in effect a res- 
onant scattering process if the small fine-structure en- 
ergy differences among the three sub-levels are ignored. 
Rather than following the photon path in detail (e.g., 
with a Monte Carlo simulation), we simply assume a sin- 
gle scattering in which the absorbed photon is re-emitted 
isotropically with the appropriate Doppler shift, and it 
either hits the star; hits the opaque, flat disk; or escapes 
the system. While this is inconsistent with the assump- 
tion of an opaque line, we find that the exact contribu- 
tion to the observed profile from the scattered photons 
has no significant bearing on our conclusions (see subse- 
quent sections), so the extra effort is unwarranted. 

The emergent spectrum at a particular viewing angle 
i is made up of photons that, either because they are 
not absorbed or they are scattered, escape into a solid 
angle bin centered on i. For a random selection of i over 
4n steradians, cosi is uniformly distributed. Considering 
five viewing angles, we choose cosi = 0.9, 0.7, 0.5, 0.3, 
and 0.1, or i = 26°, 46°, 60°, 73°, and 84°. 

4.1.1. Basic Dipolar Flow: Results 

Figure QT] shows an example of the components con- 
tributing to the emergent model profile in the case of 
R = 4i?» (0 O = 30°), and / = 0.05 (3.4 < R/R* < 4.9), 
viewed at i = 60°. The final emergent spectrum, shown 
in black, is the sum of contributions from the stellar and 
the veiling (accretion shock) continua, each shown sepa- 
rately in solid gray. The stellar contribution is the one 
with a normalized continuum level of 0.78, due to scatter- 
ing of the 4000-K continuum, while the veiling contribu- 
tion is the one with a normalized continuum level of 0.22, 
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Fig. 11. — Example scattering profile for dipolar infall with 9$ = 
30° and F = f = 0.05 (equivalently, _R = 4R* and 3.4 < R/R t < 
4.9), viewed from an angle i = 60°. The emergent profile (black 
line) is the sum of two components: the profile due to photons from 
the 4000-K stellar continuum (upper gray line) and the profile due 
to photons from the 8000-K veiling continuum (lower gray line). 
The veiling, shown in the lower left, is the ratio of the veiling 
continuum height to the stellar continuum height. Each component 
is further made up of two subprofiles: the absorption profile of the 
continuum (dashed) and the emission profile of scattered photons 
that escape toward the line of sight (dotted). 



due to scattering of the 8000-K continuum. The ratio of 
these two continua, 0.22/0.78 = 0.3, is the 1-^m veiling 
TV, also noted in the figure. Each solid gray component is 
further the sum of two subcomponents. One, shown with 
a dashed line in each case, is the absorption profile of the 
respective continuum. The other, shown with a dotted 
line in each case, is the emission profile, produced by 
scattered photons that escape toward the specified line 
of sight. The emission subcomponent is both broad and 
weak, since scattered photons can be red- or blueshiftcd 
and because for each photon absorbed, the ensuing emit- 
ted photon may hit the disk or star and not escape. Thus 
the filling-in of the red absorption by its own associated 
scattered emission is generally slight. 

Figure QTJ also illustrates an important aspect of the 
models, that the redshifted absorption in the emergent 
profile is affected differently by scattering of the stellar 
and the veiling continua. In this model, with f?o = 30° 
and i = 60°, the line of sight toward the veiling con- 
tinuum intersects the portion of the accretion flow close 
to the star where the gas velocity is high (see Fig. [T0|) . 
and scattering of the veiling continuum produces a red 
absorption that extends from 0.27 to 0.87 V BSC . In con- 
trast, the line of sight toward the stellar continuum inter- 
cepts portions of the accretion column with smaller in- 
fall speeds and a smaller velocity component is projected 
onto the line of sight. The red absorption thus produced 
ranges from ~ to 0.74 V csc and is also shallower than 
the one from the veiling continuum. The resultant ab- 
sorption profile is thus complex in shape and broader 
than either of the two individually, and in this case it 
has a maximum depth of about 20% into the summed 
continuum. 

Figures [121 [13] and [14] show the full range of model 
profiles for the three chosen values of i?o : 2, 4, and 8 i?* 
respectively (corresponding to 9q = 45°, 30°, and 20.7°). 
In each figure, the three rows show, from top to bottom, 
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Fig. 12. — Scattering profiles for dipolar infall with Ro = 2i?* (So = 45°) and F = /. Each row shows a different value of /. Within 
each row, each panel shows the profile for a different viewing angle and the corresponding l-fim veiling ry Emergent profiles (upper black 
lines) are the sum of the profiles from the stellar continuum (T = 4000 K; gray lines) and the veiling continuum (T = 8000 K; dotted lines). 
The optical veiling (at A = 5700 A) is approximately 4 times greater than ry. 



the three selections of / = 0.01, 0.05, and 0.1. In each 
row, the five panels show the profiles for the five view- 
ing angles, from i = 26 — 84°. Within each of the 15 
panels representing a unique combination of / and i, the 
final emergent profile is shown (solid black curve) along 
with the separate contributions from scattering of the 
veiling continuum (dotted curve) and the stellar contin- 
uum (solid gray curve), but not the subcomponents of 
absorption and scattered emission from each continuum 
source. 

We emphasize the following points from these three 
figures: 

1. When / = 0.01, the red absorption is dominated 
by scattering of the stellar continuum, generally showing 
small absorption equivalent widths and velocity widths. 
The profiles have a strong dependence on inclination, 
with shallow absorption at small inclinations and narrow, 
deeper, low-velocity absorption at high inclinations. 

2. The magnitude of the red absorption, measured by 
either the equivalent width or the maximum depth of ab- 
sorption, is sensitive to the parameter /. As / increases, 
there is both an increase in the veiling continuum and an 
increase in the coverage of accreting field lines projected 
in front of the stellar surface for a given Rq , enabling the 



line of sight to each point on the star to intersect more 
accreting field lines and hence yield a broader range in 
the projected velocity of the infalling gas. 

3. For a given /, the red absorption is generally 
stronger at a larger Rq (smaller 0q), since the accreting 
field lines then cover a greater range of solid angles, and 
the larger span between Ri and Rf produces a broader 
range in the gas velocity. However, inclination also plays 
a role, so that for a given / and Ro, the strongest ab- 
sorption occurs at a line of sight i that parallels the final 
part of the trajectory of the accretion flow. From the 
schematic in Figure [HJ1 it can be seen that for 6q = 45°, 
30°, and 20.7°, the corresponding viewing angle to max- 
imize the red absorption is roughly i rs 84°, 60°, and 46° 
respectively. Thus an increased 9q (smaller Rq) requires 
a higher i for a strong red absorption. This occurs be- 
cause the contribution to the absorption from scattering 
of the veiling continuum is broadest when viewed in a 
direction parallel to the flow just before it impacts the 
star. 

4. The observed emission, i.e., the part of the profile 
above the continuum, is usually weaker than the absorp- 
tion and is mostly blueshifted. Only in the extreme but 
unrealistic case when Rq = 8i?» (0q = 20.7°) with high / 
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Fig. 13. — Scattering profiles for dipolar infall with Ro = 4R* (0q = 30°) and F = /, as in Figure [T2l 



(0.05 or 0.1) and excessive Rf (13 and 35 i?» for / = 0.05 
and 0.1 respectively) does a double-peaked profile result 
when viewed close to edge on. The resulting accretion 
flow has a large solid angle, and since it is assumed to 
be in corotation with the star, the rotational broadening 
is considerable. This situation was included to complete 
our chosen parameter space and is not realistic. 

5. The veiling (r\) from the 8000-K accretion zone 
depends on /, 6q and i. The dependence on / is ob- 
vious, since r\ scales almost linearly with / (eq. [4]). 
The dependences on 6$ an< i * arise through their influ- 
ences on the projected area of each continuum source. 
The parameter 9q signifies the orientation of the veil- 
ing continuum, hence its direct effect on the veiled area 
viewed. Although less sensitive to the viewing angle, the 
projected stellar continuum area changes because of the 
presence of the disk extending from Ri (dependent on 8q 
and /) to infinity. For both 9 = 20.7° and 30°, r Y at 
a given / drops monotonically as i increases from pole- 
on to edge-on, by a factor of 5 and 3 respectively. At 
#o = 45°, ry varies less, dropping by a factor of 1.6 from 
i = 26° to 73°, then increasing slightly toward i = 84°. 
For example, when R Q = 4R* (0 O = 30°) and / = 0.1, 
ry ranges from 0.4 to 1.15 and ry ranges from 1.5 to 4.5 
with viewing angle. 

In sum, a strong red absorption extending to high 



velocities, like those observed, requires / (and thus r\) 
to be large so the contribution from scattering of the 
veiling continuum is enhanced and the angular extent 
of the flow on the star is increased. Strong, broad 
absorption is also more likely when Rq is large and the 
line of sight parallels the accretion flow close to the star. 
The relation between absorption magnitude and r\ is a 
crucial test of the dipolar accretion model, as we will 
show in the following subsection when we compare the 
observations to our model profiles. 



4.1.2. Basic Dipolar Flow: Comparison to Observations 

In comparing our models with observed profiles, we fo- 
cus on the red absorptions, since the small emission at 
blueward velocities expected from scattering in the fun- 
nel flow will often be overwhelmed by additional sources 
of emission, such as scattering and in-situ emission from 
a wind. The red absorptions are evaluated in context 
with the observed veiling, which is the basis for evalu- 
ating /. It is immediately apparent that there is a mis- 
match between the model profiles in Figures [T^l EH and 
[T4l and the observed spectra in Figure O since the ma- 
jority of CTTS are known to have / < 0.01 (CG) while 
model sequences for / = 0.01 (implied ry ~ 0.06 and 
ry ~ 0.25) have shallow and/or narrow red absorptions 
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Fig. 14. — Scattering profiles for dipolar infall with i?o = &R* (0o = 20.7°) and F = f, as in Figure |j"2l 



bearing little resemblance to the ensemble of broad and 
deep observed He I A10830 profiles. 

A more explicit demonstration of the limitation of the 
models can be made from a quantitative comparison be- 
tween the equivalent width and veiling for model and 
observed profiles. This comparison requires normaliz- 
ing the observed profiles to their respective escape ve- 
locities, since in the models all velocities are in units of 
the escape velocity. The normalized equivalent width, 
W' x = W\/V Qac , has an intuitive interpretation: it is sim- 
ply the fraction of the continuum absorbed between rest 
and the escape velocity, with a value of 1 indicating to- 
tal absorption over the entire range. Figure [T51 compares 
the normalized red absorption equivalent width W' x to ry 
and ry for both models and observations of the reference 
sample. In the models we have assumed that the excesses 
at both Y and V arise from an accretion shock that emits 
an 8000-K blackbody continuum. This is known to be a 
valid assumption for optical veilings, and while the ry 
data points are not simultaneous with the observed ab- 
sorption profiles, the fact that all but one of the objects 
with ry = also have low ry indicates that this is a rea- 
sonable approach. Unless the ry values for these objects 
were all a factor of 5 to 10 higher when the He I A10830 
profiles were obtained than when the HEG data were ob- 
tained, the two panels together clearly indicate that only 



a fraction of the observed data lie within the realm of 
the model results: those with weak red absorption and 
small veiling or those with modest red absorption and 
intermediate veiling. There is a glaring discrepancy be- 
tween models and observations for stars with large W' x 
and small ry. 

The values for W' x and adopted escape velocities are 
listed in Table [5] along with two additional properties 
of the red absorption: the normalized width and the 
depth. The normalized full-width at quarter-minimum, 
FWQM', is the width measured at one quarter of the 
absorption minimum as a fraction of the escape veloc- 
ity. The depth of the absorption component at 0.75 V esc 
normalized to 100%, Do. 75, was chosen since it is sen- 
sitive to the infall geometry close to the star. As with 
W x , a number of stars indicate a discrepancy with the 
basic dipole model, where objects with small veiling fre- 
quently have both FWQM' and Z?o.75 much larger than 
can be accounted for with the models. This is illustrated 
in Figure [TBI which shows the comparison of the velocity- 
normalized equivalent width, velocity-normalized line 
width, and high-velocity depth to 1-fim veiling for ob- 
servations and models. 

This comparison demonstrates that a fraction of CTTS 
with subcontinuum red absorption at He I A10830 have 
red absorptions too strong to be accounted for by mag- 
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TABLE 5 

Measurements of Red Absorption in 

Velocity-Normalized He I A10830 Profiles 
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<r u > 



Fig. 15. — Comparison of the red absorption equivalent width 
(normalized to the escape velocity) to the 1-fim veiling (left) and 
the non-simultaneous average optical veiling (right; from HEG) for 
basic dipolar models and the profiles from the reference sample. 
The model properties appear as lines connected by symbols. Each 
symbol type is for a different Ro / #0i with circles for Ro = 2R*, 
asterisks for Ro = 4iJ*, and diamonds for Ro = 8-R*. Each line 
type is for a different viewing angle, with solid black for 26° , solid 
gray for 46° , dotted for 60° , dashed gray for 73° , and dashed black 
for 84°. Along a line, symbols indicate filling factors / = 0.01, 
0.05, and 0.10, always increasing toward increasing veiling. Since 
the veiling axes are logarithmic, stars with no detected l-/im veiling 
are placed at ry = 0.025. 
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Fig. 16. — Comparison of 3 properties of the red absorption to 
the 1-fj.ra veiling for basic dipolar models and the 21 profiles from 
the reference sample. In each column, the observed parameters 
are the same, but model parameters appear only for the indicated 
Ro I 8o combination. The model parameters for the three filling 
factors / = 0.01, 0.05, and 0.10 at a particular viewing angle are 
connected by lines, and the correspondence between line type and 
viewing angle is the same as in Figure [TBI As before, stars with no 
detected 1-fj.m veiling are placed at ry = 0.025. 
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0.18 
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UZ Tau W 
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0.54 


3 


V836 Tau. 
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0.11 


0.39 
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YY Ori . . . 


[430] 


0.13 


0.49 


11 



Note. — For the reference sample only. Col. 2: Adopted 
escape velocity (km s _1 ); Col. 3: Equivalent width of 
velocity-normalized absorption (dimensionless); Col. 4: 
Full-width at quarter-minimum of velocity-normalized ab- 
sorption (dimensionless); Col. 5: Depth at 75% of the es- 
cape velocity as a percentage of the continuum. 
a Brackets indicate V ro d/V eS c > 1, so we assume V csc » = 
V rcd /0.9. 

netospheric accretion in a basic dipole, where the filling 
factor of the flow on the stellar surface F is equivalent 
to the filling factor of shocked gas at the terminus of ac- 
creting field lines /. The observations that present the 
greatest challenge to the model are those in which the red 
absorption is strong {W' x > 0.1) but the veiling is weak 
(ry < 0.1). This conclusion is robust, since the models 
have been constructed to produce maximal red absorp- 
tion for a given Rq and /, in that the A10830 transition is 
assumed to be optically thick and the thermal/turbulent 
broadening has a generous 10 km s _1 half- width. 

This conclusion is not compromised by the choice of 
8000 K for the temperature of the shock-heated pho- 
tosphere. This value corresponds to the low end of 
the temperature range derived fro m modeling the SEDs 
of observed continuum excesses (Harti gan et alj I1991L 
iGullbring et~aTlll998t I Johns-Krull et al.ll2000O f higher 
temperatures were assumed, the veiling for a given / 
would be even larger, worsening the agreement between 
the models and the observations. If we adopted the low- 
est temperature allowed by the SED models of the optical 
continuum excess, T ~ 6000 K, the associated veiling for 
a given / would be reduced by ~ 2 at Y and ~ 3 at 
V. Figure [15] demonstrates that shifting all the model 
results to the left by a factor of 2 in ry or 3 in ry is still 
insufficient to account for the strong absorptions and low 
veilings. We thus conclude that those profiles with strong 
absorptions and small veilings lie outside the realm of 
model results for self-consistent dipole flows. 

4.2. Dilution 

A simple way to keep the veiling small and yet have 
the accretion flow project a broad velocity range in front 
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of the star is to let the flow arise over a large range of R 
(thus impacting the star over a large range of 9) but to 
fill the whole enclosed volume only dilutely with accret- 
ing gas. We now distinguish between F, the fractional 
surface area on the star over which the magnetosphcric 
footpoints are distributed, and /, the fractional surface 
area on the star occupied by accretion shocks at the base 
of field lines that carry accreting gas. We define /' = f/F 
as the fraction of F occupied by all the accretion shocks. 
With enough dilution, i.e., /' sufficiently small, F can be 
large enough to provide the areal coverage over a large 
velocity range that is necessary for a broad and deep red 
absorption, while / = Ff can remain small, as required 
to produce a low veiling. One way to achieve this is 
to postulate a large number of narrow accretion stream- 
lets spatially separated from one another that together 
impact only a fraction /' of the outlined area F. (We as- 
sume the many accretion shocks are dispersed randomly 
throughout F.) Then, with an intrinsic thermal or tur- 
bulent line broadening of ~ 10 km s _1 associated with 
each streamlet, photons from the star can intersect a 
sufficient number of streamlets such that the continuum 
(stellar or veiling) will be absorbed over the full velocity 
range specified by the parameter F as though the whole 
volume were filled. 

The concept of many accretion streamlets dilutely 
filling a volume has the additional advantage of offer- 
ing a credible explanation for how the lower level of 
A10830 (2s 3 S) is populated over all streamlines. With 
the difficulty of maintaining a temperature high enough 
(> 2 x 10 4 K) for collisional excitation to the 2s 3 S level 
in a freely falling gas, it is likely that photoionization is 
the excitation mechanism. Then, if the source of ioniz- 
ing photons is the accretion shock itself, the much smaller 
shocked area of an individual streamlet within a diluted 
flow, as compared to the shocked area of a single undi- 
luted flow, will enable more ionizing photons to escape 
from the sides and ionize the gas in other streamlets, 
even at positions far from the star. Or, if the domi- 
nant source of ionizing radiation is located away from 
the streamlets (e.g., the stellar corona), these photons 
will be able to penetrate into the volume and ionize in- 
dividual streamlets as opposed to ionizing just the skin 
of a single completely filled accretion flow. Thus, many 
narrow streamlets dilutely filling a large volume not only 
yield a deep red absorption from the large coverage area 
over a broad velocity range of infalling gas, but they also 
readily account for the ionization of gas at each location 
in the flow to produce an optically thick A10830 transi- 
tion over the whole velocity range. 

4.2.1. Profiles for Wide, Dilutely Filled Flows 

We compute scattering profiles for diluted dipole flows 
for the same 3 geometries shown earlier, with (Rq, 0q) 
pairs of (2 i?», 45°), (4 #*, 30°), and (8 #*, 20.7°). We 
introduce a wider range in F, from 0.01 to 0.2, although 
now all models have / — 0.01, corresponding to a range 
in /' from 1 to 0.05. The resulting profiles are shown 
in Figure \T7\ and the model parameters are listed in the 
lower portion of Table |H In the figure, the 3 columns 
correspond to the 3 Rq values and each row is a common 
value of /'. The degree of dilution increases downward in 
the figure, with the case for no dilution shown in the top 
row (/' = 1 and f — F) repeated from Figures [121 [T3l 
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Fig. 17. — Red side of scattering profiles for a series of "diluted" 
dipoles all with / = 0.01 (ry < 0.11) but with F ranging from 
0.01 (top row) to 0.2 (bottom row). From left to right, columns 
correspond to Rq = 2, 4, and 8 i?*, and the range of R for each 
F is specified. The profile sequences for each panel correspond to 
viewing angles of 26°, 46°, 60° (black), 73°, and 84°. The top row 
with F = f = 0.01 corresponds to the models shown in the upper 
rows of Figures |12I 1131 and 1141 For comparison with observations, 
dotted horizontal lines mark depths of 20% and 30% and crosses 
mark a depth at V/Vme = 0.75 (A).7s) of 10%. 

and 1141 In subsequent rows the dilution grows as F in- 
creases to 0.05, 0.1, and finally 0.2. Each panel shows the 
superposed profiles for all 5 viewing angles for each i?o, 
/' (or F) combination and the corresponding inner and 
outer radii of the accreting volume, Ri < R < Rf. Since 
the effect of the viewing angle on the profile morphology 
is roughly independent of dilution, the individual viewing 
angles can be identified by referring to the earlier figures. 
We highlight the i = 60° profile with a darker line, since 
this is the most probable viewing angle. 

Although all profiles in Figure Q~7] are for / = 0.01, the 
associated veilings differ slightly because the accretion- 
heated area is distributed differently for different values 
of F and 8 Q , leading to slightly different projected 
areas. Nonetheless, in all cases, ry < 0.11. For these 
small veilings, the absorption, while quite strong when 
there is significant dilution, is almost entirely due to 
scattering of the stellar continuum, in contrast to the 
undiluted models where large veilings and scattering of 
the veiling continuum were necessary to produce strong 
absorption. As dilution increases for a given i?o, the red 
absorptions become increasingly strong and broad, due 
to the increased areal coverage over a broader range of 
velocities as the interval between Ri and Rf increases. 
For example, the maximum penetration depth of the red 
absorption into the continuum, -D maX 7 increases from 
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10% to 30% for Ro = 2i?* between an undiluted and a 
/' = 0.05 flow. As before, larger Rq also increases the 
areal coverage and thus the depth and breadth of the 
absorption: For Rq = 4i?» and F — 0.2, D max reaches 
to 50% of the stellar continuum for all viewing angles. 
In the unrealistic case of Rq — 8i?» and F = 0.1, where 
Rf extends to 35 i?*, D max can be 70% of the stellar 
continuum. However, for flows confined to maximum 
sizes on the order of corotation, the deepest penetrations 
are about 50% of the stellar continuum. 



4.2.2. Further Comparison to Observations 

By introducing the concept of a diluted dipole, where 
field lines carrying accreting gas only dilutely fill the vol- 
ume occupied by a wide magnetosphere, we can simulta- 
neously generate deep and broad red absorption features 
while maintaining small filling factors for hot accretion 
shocks with / ~ 0.01. This is the empirical regime in 
Figure [TBI where dipolar models with f = F were unable 
to account for stars with both strong absorption (large 
W' x ) and very low veiling (ry ~ 0). We compare observed 
and model profiles for a few individual stars in Figure fTSl 
in 4 cases for undiluted, fairly narrow dipoles for stars 
with ry ranging from to 0.4 and in 2 cases for wide, 
diluted dipoles with ry ~ 0, where we have rescaled the 
model profiles to the escape velocity of each star. Since 
we have not computed a large grid of models, the mag- 
netospheric properties listed for each fit are not intended 
to be predictions for a particular star. However, this fit- 
ting procedure shows that weaker red absorptions can be 
reasonably described by basic undiluted models with a 
small range of origination radii in the disk, where veilings 
ry from to 0.4 can be consistently modeled with an ap- 
propriate choice of /, and the red absorption can include 
scattering contributions from both stellar and accretion 
shock continua. Similarly, strong red absorptions in stars 
with low veilings can be well fit by dilutely filled flows 
with small / but a wide span of origination radii in the 
disk, resulting in a large projected area of accreting gas 
for the scattering of the stellar continuum. 

The overall applicability of the diluted dipolar model 
can be appreciated by comparing the model profiles from 
Figure IPTl to the ensemble of observed helium profiles for 
those stars with ry < 0.1 and thus / ~ 0.01, where the 
effect of scattering from a hot accretion shock will be 
inconsequential and the properties of the red absorption 
will be shaped almost entirely by scattering of the stellar 
continuum. To effect this comparison in a general way, 
rather than focusing on individual stars, in Figure [151 we 
plot superposed observed profiles for the redward side of 
He I A10830, each normalized to their respective escape 
velocity and separated into 3 groups on the basis of their 
depths both at 0.75 V esc (Do. 75) and at maximum absorp- 
tion (D max )- To aid in the comparison, both Figures IT71 
and [T9l denote depths for Do 75 = 10% and D max = 20% 
and 30%. 

In Figure [19] the left panel contains the 3 shallowest 
profiles, with D max < 20% and D0.75 < 10%. Compared 
to the predicted profiles in Figure \TT\ the model flows 
that most resemble such broad but shallow profiles have 
Ro ~ 2R* and F < 0.1, although some viewing angles 
for larger flows with relatively small areal coverage of 
magnetic footpoints, F < 0.05, could also apply. Our 
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Fig. 18. — Examples of least-squares fits of dipolar model pro- 
files (gray) to selected observations (black), with the stellar escape 
velocities marked by short vertical lines. The top four panels use 
dipoles with F = /, and the corresponding /, ry, and i are shown. 
The bottom two panels have strong red absorptions and no de- 
tected veiling; they are well fit by extended dilute dipolar models 
(Ro = 4fl*, F = 0.2 for GI Tau; Ro = 4R», F = 0.1 for V836 
Tau). Since processes other than scattering by the accretion flow 
can be important at low velocities, points with V/V esc < 0.1 are 
ignored in the fitting procedure. 
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Fig. 19. — Superposed He I A 10830 lines from the reference 
sample for the 13 stars with ry < 0.1, appropriate for modeling 
with dilute dipole flows with / = 0.01. Only the red half of the 
profile is shown, normalized to the individual escape velocity of 
each star. Profiles are grouped by -Do. 75, the penetration depth 
into the continuum at V/V CB c = 0.75, and by -D ma x, the maximum 
penetration into the continuum. For comparison with models, the 
dotted cross in each panel marks Do. 75 = 10%, and the dotted 
horizontal lines mark depths of 20% and 30%. 

coverage of parameter space is not exhaustive, but it is 
clear that for the broad but shallow red absorptions the 
range of radii over which the accretion flow leaves the disk 
is narrow, corresponding to a fairly small area on the 
star for the magnetospheric footpoints, but still larger 
than 1%. The central and right panels contain the 10 
deeper profiles among the stars with low veiling, where 
Dmax ranges from 30 to 60%. Model flows that produce 
deeper profiles generally have significant areal coverage 
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of magnetic footpoints F, as seen in Figure [T7l where the 
accretion flow leaves the disk over a wide range of radii, 
impacting the star over a wide range of angles, in some 
cases with magnetic footpoint coverage up to 20% of the 
stellar surface area. This is considerably larger than has 
been modeled in previous work on magnetospheric infall. 

Even with wide diluted flows, the profiles of the 3 stars 
in the right panel of Figure [19] (AA Tau, DK Tau, DN 
Tau) are a challenge to explain under the constraints 
of a dipolar geometry. These profiles not only have 
-D max > 30% but also have -Do. 75 > 10%, with the caveat 
that errors in escape velocity may be up to 20%. From 
the models explored in Figure I17[ flows with very wide 
extents, leaving the disk over a range of radii from a few 
i?» to beyond corotation and viewed fairly close to pole- 
on, are required to produce profiles with -Do. 75 > 10%. 
Rather than postulate an enormous dipolar flow with a 
polar viewing angle (which is clearly not the case for, at 
least, the edge-on source AA Tau), in the next section 
we will explore an example of a non-dipolar geometry to 
find a more plausible explanation for these three obser- 
vations. 

Only 8 CTTS in the reference sample have ry > 0.1, 
such that the properties of the He I A10830 red absorp- 
tion may be affected by scattering of continuum photons 
from the hot accretion shock. One of these is DR Tau, 
where the high ry = 2 implies / « 0.24 (eq. [4]), which 
as shown in Section 4.1 would yield a red absorption at 
least an order of magnitude stronger than the observed 
W' x — 0.05. As will be addressed in Section 5, we suspect 
that in this case the red absorption has been filled in by 
a wind exterior to the accretion flow. 

4.3. Diluted Radial Flows 

We have identified the 3 stars in the right panel of 
Figure[H AA Tau, DK Tau, and DN Tau, as difficult to 
explain with scattering in a dipolar geometry due to their 
absorption depths at velocities in excess of 0.5 V csc . In 
a dipolar flow, the impact velocity at the stellar surface 
depends on the polar angle 9, which is determined by 
the initial distance of infall R (eq. [3]), such that the 
impact velocity is greatest when 9 is near the pole (i.e., 
R is large) and diminishes as 9 approaches the equator 
(i.e., R becomes small). Thus if 9 is small enough, high 
impact velocities will result, although flows with small 
9 become highly curved and pinched as they reach the 
star (Fig. [i"0|) , resulting in small areal coverage and thus 
a shallow absorption profile at the highest velocities. 

We investigate radial infall trajectories as an alter- 
native geometry that could produce deep absorption at 
high velocities. In all aspects except the geometry, ra- 
dial models have the same assumptions as our dipolar 
models except that we have not included rotation. The 
axisymmetric flow begins at some distance from the star 
i? max , and it falls radially toward the star, impacting 
the stellar surface between polar angles 9\ and #2 in one 
hemisphere and between 7r — 9\ and 7r — 9i in the other. 
The fractional surface of the star spanned by the accre- 
tion flow, F, is cos ^i — cos 02, and the shocks within 
this region together occupy a fraction /' of the area F, 
so that / = Ff. The disk truncation radius is a free 
parameter, but we set it equal to i? m ax, which is 8 -R* 
in all radial models. Figure [501 shows scattering profiles 
from two radial geometries at 5 viewing angles. In the 



top row, the impact region extends from 9\ = 66.4° to 
9i = 78.5°, while in the bottom row, the impact region 
extends from 9\ — 78.5° to the equator. In both cases, 
F = 0.2 and / = 0.01. As expected, the absorption is 
strongest for a viewing angle within the confines of the 
flow (i.e., 9\ < i < 62), and the profile becomes a nearly 
symmetric emission profile (assuming axisymmetry and 
no rotation) for views close to pole-on. When the view- 
ing angle is aligned or nearly aligned with the column 
of absorbing gas, each radial model can produce the ob- 
served range of absorption depths at high velocities, with 
-Do. 75 > 10% for profiles with i > 60° in the top row and 
i > 73° in the bottom row. 

We are not advocating radial infall starting from a 
large distance, and thus, the profile sequences in Fig- 
ure [20] are not expected to be realistic for the whole 
velocity range. However, the requisite deep absorption 
at high velocities, resulting from material moving faster 
than ~ 2/3 V CS c, all arises inside of about 2 i?». Thus 
the message from these calculations is that the accretion 
stream only need move in a radial trajectory, i.e., become 
less curved than a dipole, as it nears the star. 

The effectiveness of radial infall trajectories for ma- 
terial near the star in accounting for the high-velocity 
absorption in AA Tau, DK Tau, and DN Tau is shown 
in Figure [2X1 The figure shows model and observed pro- 
files where (1) profiles are inverted so the vertical axis is a 
measure of the minimum stellar coverage fraction at each 
velocity and (2) only velocities in excess of 0.5 V csc are 
plotted. The regime of diluted dipolar models with the 
largest -Do. 75 is shown with dark and hatched shading, 
while the regime of flows with radial trajectories for gas 
near the star is shown with light shading. The dark shad- 
ing is for the best case from our diluted dipolar models 
for a flow contained entirely within the corotation ra- 
dius: F = 0.1 originating between 3.0 and 6.2 i?„ in the 
disk and viewed from i — 46°. Although 8 R* is a more 
typical corotation radius, extending the flow out to this 
distance would not produce much additional absorption. 
The hatched region is for a diluted dipole that allows 
field lines extending out to nearly twice the corotation 
radius to participate in the flow, where the dashed lines 
are for the case R = 4R*, F = 0.20 (2.4 < R/R* < 15), 
seen from two viewing angles, i = 26° and i — 46°. Al- 
though the latter two extreme dipolar models come close 
to producing sufficient absorption at high velocities, sig- 
nificant accretion beyond corotation is likely not phys- 
ical. In contrast to the dipole trajectories, the regime 
of the four radial models from Figure [20] with viewing 
angles nearly aligned with the infalling gas easily con- 
tains the observed stellar coverage fraction from 0.6 to 
0.85 V csc , with no requirement that the flow originate at 
radii beyond corotation. The realistic situation is likely 
to involve some complex magnetic field topologies with 
trajectories approaching radial as they near the star. 

5. DISCUSSION 

5.1. Implications of Diluted Funnel Flows 

The high opacity and resonance scattering properties 
of He I A10830 enable the geometry of magnetospheric 
accretion to be probed via absorption of gas seen in pro- 
jection against the star, in contrast to previous studies 
that rely on the morphology of emission lines. Under 
the assumptions that the flow is an azimuthally symmet- 
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Fig. 20. — Scattering profiles for diluted radial infall in non-rotating, azimuthally symmetric flows that begin at 8 R* and impact the 
star over a range of polar angles 8 that encompasses F = 20% of the stellar surface area. Accreting field lines and their accretion shocks fill 
only 5% of F (i.e., / = 1%), with 1-^tm veilings ry as listed. In the top row, the flow impacts the star over the range 66° < 8 < 78°, while 
in the bottom row, the flow impacts the star over the range 78° < 9 < 90°. The same five viewing angles are used as in previous figures. 

since the strongest and broadest He I A 10830 absorp- 
tions are seen in stars with little or no l-/im veiling, 
these red absorptions must instead arise almost solely 
by scattering of photospheric radiation. Achieving the 
observed breadth and depth of the absorption requires 
a large angular coverage of the stellar continuum in the 
azimuthal direction over a wide range of velocities for 
many stars, with areal coverage in footpoints on the star 
of F = 10 — 20%. We suggest that the required com- 
bination of wide flows and low filling factors of hot gas 
is a result of accretion in many narrow streamlets, each 
of which may have a dipolar configuration but which to- 
gether only fill a small fraction of the enclosed volume. 
We have explored the case where the streamlets are uni- 
formly distributed through the accreting volume, pro- 
ducing wide, dilutely filled flows that reconcile the need 
for absorption over a broad range of velocities with filling 
factors of hot gas / < 1%, as observed (CG). 

Earlier studies also imply a discrepancy between the 
areal coverage F of magnetospheric footpoints and the 
filling factor of hot accretion shocks /. For example, 
magnetospheres with f — F — 8% were invoked to 
model hydrogen lines arising from accretion flows in 
order to pro duce sufficient line flux es and mass accre - 
tion rates ([Symington et al1 l2005: K urosawa et al J 2006) . 
The seminal sequence of pape rs modeling hydrogen l ine 
formation in funnel flow s from lHartmann et al.l (|1994f ) to 
iMuzerolle et al.l (|2001[ ) also required filling factors that 
were larger than predicted by SED modeling of contin- 
uum excesses to account for observed emission line lu- 
minosities. The notion of accretion via streamlets that 
dilutely fill a large volume is a straightforward way to 
reconcile this discrepancy, simultaneously allowing large 
field sizes and small shock filling factors. 




Fig. 21. — High- velocity tails of observed and model profiles for 3 
stars with the largest values of -Do. 75 an d ry = 0, inverted to show 
the minimum fraction of the star occulted by infalling material at 
each velocity. Dark shading indicates the regime of diluted dipolar 
models with / = 0.01 and infall contained entirely within a typical 
corotation radius, marked by the profile (dotted line) with F = 0.1, 
3.0 < R/R* < 6.2, and i = 46°. The hatched region indicates the 
extension when dipolar field lines out to ~ twice the corotation 
radius participate in infall, marked by the profiles (dashed lines) 
with F = 0.2, 2.4 < R/R* < 15, and i = 26° or 46°. Light shading 
shows the regime for profiles formed in diluted radial infall with 
F = 0.2 and / = 0.01. For both radial geometries in Figure [20l 
two profiles (solid lines) with i close to the infall angle are shown. 

ric dipole and helium is sufficiently optically thick that 
all incident l-//m radiation is scattered, we have illus- 
trated the sensitivity of the red absorption to both the 
angular extent of the magnetosphere and the filling fac- 
tor of hot gas from the accretion shock /. If / exceeds a 
few percent, the hot spot will be an important contrib- 
utor to the scattering of the l-/im continuum; however, 
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Although our model invokes diluted accretion flows in 
widely distributed streamlets of gas, an alternate sce- 
nario for diluted acc retion is the one sugge sted by the 
MHD simulations of iRomanova et all (|2004f ). where in- 
ternal structure within the accretion flow gives a mass 
accretion rate (and a corresponding blackbody contin- 
uum temperature) that is highest at the interior and falls 
off toward the sides. Although this scenario can also, 
for a large F, produce a smaller veiling from the area- 
weighted blackbody continua than the undiluted / = F 
case, the advantage to widely dispersed streamlets is that 
they provide a facile means for ionizing radiation to pen- 
etrate to most of the infalling gas, since distributed ac- 
cretion shocks with small individual areas would allow 
ionizing photons produced in each shock to escape more 
easily from the sides and ionize helium at other loca- 
tions. Another consequence of such distributed accretion 
shocks is that photons from the shocks emitted toward 
the star would be incident on a larger area of the pho- 
tosphere than for a single shock with the same /. This 
may invalidate the usual assumption of a plane-parallel 
geometry for the radiative transfer of photons with the 
effect that, independent of the internal structure within 
an individual streamlet, the resultant veiling continuum 
would encompass a range in blackbody temperatures. In 
a wide flow where dilution is somewhat uniform, there 
will be many separate shocks with a range of blackbody 
temperatures surrounding them. There may be some 
observational support for this phenom enon in that the 
veiling continuum longward of 0.5 /im ( Basri fc Batalhal 
ll990tJWhite k. Hillenbra nd 2004; EFHK) is broader than 
the single 8000-K blackbody that is a good match to the 
excess at shorter wavelengths (CG). 

Constraints on the angular extent of accreting gas and 
the location in the disk where infall originates are rele- 
vant to models for disk locking and wind launching. Al- 
though there are some cases where He I A 10830 profiles 
resemble those expected from viewing an accretion funnel 
restricted to a narrow origination around the corotation 
radius, the suite of profiles expected from viewing this 
magnetic topology from all inclination angles is not con- 
sistent with the observations. The most extreme deep 
and broad absorptions instead require infall spanning a 
wide extent of origination radii, from a few i?* out to 
at least typical corotation radii of 6 to 8 i?* if the flows 
are dipolar. For other magnetic field configurations, such 
as a tilted dipole or a multipolc field, significant red ab- 
sorption need not require such a wide range in initial 
infall distances. In general, the depth of the red absorp- 
tion is governed more by the range of impact latitudes 
than by the range of initial radii; only for a dipolar flow 
aligned with the rotational axis are the two ranges so 
closely linked. For example, in an aligned dipolar flow 
with Rq = 4R* and F = 0.1, the range in impact latitude 
of 23.6-35.3° corresponds to a range in initial radius of 
3.0-6.2 i?». In a more complex magnetic configuration, a 
comparably wide range of impact angles could produce a 
strong red absorption without the need for such a large 
range of initial radii. 

The necessity for a dilutely filled flow does imply that 
there is not a sharp delineation on the disk for accretion 
onto the star. It likely indicates a very inhomogeneous 
field structure at large distances, with many local pockets 
distributed over a broad radial range on the disk giving 



rise to accretion streamlets. Since our analysis assumes 
axial symmetry in a set of nested dipolar flows, the con- 
straints that the breadth and depth of the red absorption 
place on the angular extent of the accreting gas are even 
more extreme if, as likely, accretion channels arc in re- 
stricted azimuth zones. Furthermore, there are some red 
absorptions that are so deep at velocities > 0.5 V osc that 
a dipole morphology is inadequate, even when arising 
from 2 i?* to the corotation radius. In these cases we 
find that radially directed infall can achieve the requisite 
depth of absorption, although other topologies that re- 
sult in a large covering factor of the star at the highest 
velocities can likely be constructed. 

Recent Doppler tomographic maps of the CTTS V2129 
Oph and BP Tau, based o n circular polarization of Ca II 
(|Donati et al.l l2007. 2008|), reveal the locations of accre- 
tion hot spots on the stars. The spots span quite a 
broad latitude range (extending roughly from the pole 
to 45°) but a very narrow azimuthal range. The narrow 
azimuthal range implies that the detection of He I A10830 
red absorption requires an opportune time at which the 
accretion spots are directly in view. This situation is 
consistent with the result that BP Tau, a mildly accret- 
ing CTTS included in our He I A10830 survey, did not 
show any He I A10830 red absorption on the two occa- 
sions we observed it. At present, there are not enough 
tomographic data to see how consistent this pattern is 
among a range of accreting stars, although our detection 
of subcontinuum red absorption in 21 of 38 CTTS, in- 
cluding 20 of 29 stars (and 37 of 56 total spectra) with 
ry < 0.5 (ry < 2), would imply a large azimuthal cover- 
age by the accretion spots. However, we note that even 
in the two stars with tomographic maps, there is a possi- 
bility that accretion impacts the stars over a wide range 
of longitudes. Donati et al. attribute only 2/3 of the 
He I A5876 emission but all of the He I A5876 circular po- 
larization to accretion spots, based on the fraction of the 
emission that shows rotational modulation compared to 
that which is time-independent. The time-independent 
component, responsible for 1/3 of the He I A5876 emis- 
sion, is attributed to a chromosphcric component dis- 
tributed uniformly over the stellar surface. However, 
since non-accreting WTTS show either very weak or, 
more commonly, no He I A5876 emission (jBeristain et al.1 
2001), it would appear that TTS chromospheres are not 
significant contributors to this line. Instead, the time- 
independent component may be from more widely dis- 
tributed accretion shocks that cover a broader range of 
longitudes. 

5.2. Absence of He I X10830 Red Absorption 

In this paper we have focused on the 21/38 CTTS that 
show redshifted absorption in He I A10830 at least once 
in an observational program with sporadic time cover- 
age. Clearly the absence of He I A10830 red absorp- 
tion is also important in constraining the topology of 
magnetospheric accretion. An important point is that 
He I A10830 red absorption is rarely seen among CTTS 
with the highest l-/im veiling (1/25 observations; see Sec- 
tion 3.1). Of the 9 stars in the EFHK survey in this 
category, the only one that showed redshifted absorp- 
tion, on 1 of 4 occasions, is DR Tau. We suspect that 
in all 9 of these stars, emission from a wind exterior to 
the accretion flow, instead of from the flow itself, is fill- 
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ing in any redshifted absorption that may be present. 
If in-situ emission from the funnel flow were significant, 
it would be difficult for it to fill up the absorption at 
the red edge of the profile, since the geometry of the 
funnel flow results in smaller volumes at higher veloc- 
ities, producing centrally peaked emission profiles that 
fall off rapidly toward both blue and red high velocities 
(sec the contribution to the emission from scattering of 
the stellar continuum in Fig. ITT)) . The near absence of 
red absorption among these stars instead calls for a sit- 
uation in which the redshifted absorption, if present, is 
filled in completely. In the case of DR Tau, it is clear 
that weak red absorption, confined to high velocities, is 
visible when the emission from the P Cygni wind profile 
is weakest (see Fig. [7]). Among the other stars in this 
high- veiling group, all have either broad blue helium ab- 
sorptions indicative of viewing through a stellar wind or 
strong helium emission interpreted to arise in a conical 
stellar wind viewed obliquely (see KEF and EFHK). Ei- 
ther of these contributions to redward emission would be 
sufficient to fill up even a strong red absorption, provided 
the wind is optically thick and exterior to the accretion 
flow. 

When profiles from both He I A10830 and He I A5876 
are considered, the evidence suggests that He I A10830 
red absorption is rare or absent in CTTS with large 1- 
yum veiling not primarily because the absorption is be- 
ing filled in by wind emission but more because the ge- 
ometry of the funnel flows is altered compared to that 
of low- veiling CTTS. This inference is drawn from a 
study of He I A5876 pr o files a nd optical veiling pre- 
sented in iBeristain et al.l (|200lD , which includes many 
stars in common with EFHK. They found that CTTS 
whose He I A5876 profiles showed only a narrow compo- 
nent, consistent with formation in post-shock gas from 
an accretion shock, show an excellent correlation be- 
tween the strength of narrow-component helium emis- 
sion and optical veiling. In contrast, CTTS whose 
He I A5876 profiles show a contribution from a broad 
component have reduced or absent emission from a nar- 
row component relative to stars of similar optical veiling. 
While it might appear otherwise, this is not an esoteric 
point! It suggests that CTTS with strong stellar winds 
and high optical and l-/im veilings may have crunched 
or otherwise-altered magnetospheres resulting in weak 
narrow-component emission from a hot accretion shock, 
and there is a significant contribution to the veiling con- 
tinuum from another source. In contrast, CTTS with- 
out strong stellar winds (many of which show disk wind 
profiles at He I A10830; see KEF) have extensive mag- 
netospheres carrying accreting gas to the star, and hot 
accretion shocks are the dominant contributor to their 
optical veiling. We anticipate being able to test this 
suggestion shortly, following analysis of simultaneously 
obtained spectra extending from 0.4 to 2.2 /im. 

A second point regarding the frequency of He I A10830 
red absorption is that, in contrast to CTTS with high 
l-/xm veiling, red absorption is commonly seen in stars 
with lower veiling (37 out of 56 spectra for ry < 0.5; 
see Section 3.1). Among this group, some objects (e.g., 
TW Hya and CY Tau; see Fig. [7]) clearly show reduction 
of the red absorption as the emission, likely that of a 
stellar wind as indicated by the strong P Cygni profile, 
increases. In such stars the appearance and disappear- 



ance of the red absorption is likely due, at least in part, 
to filling in by an exterior stellar wind, as in DR Tau. In 
others (e.g., V836 Tau and GK Tau; Fig. [7]), the weaker 
helium emission could arise simply from scattering in the 
funnel flow, and the absence of red absorption may indi- 
cate viewing at an azimuth with no funnel flow activity. 
Azimuthal asymmetry in the funnel flow is also the likely 
explanation for the strongly variable red absorption mor- 
phology in objects such as AA Tau and DK Tau (Fig. [7]). 
The possibility that red absorption may be partially 
filled in, cither by in-situ emission from the accretion 
flow or by scattered or in-situ emission from a wind, im- 
plies that the true magnitudes of some red absorptions 
are stronger and their constraints on the flow structure 
stiffer than the observations indicate. Further, since red 
absorption can be completely filled in by in-situ wind 
emission or, in some cases, not be observed at all due to 
azimuthal inhomogeneities, He I A10830 red absorption 
is likely more pervasive among CTTS than is already 
apparent. 

5.3. Size and Structure of the Accretion Flow 

Inferences to date on the physical extent of accretion 
flows have largely relied on models positing that hydro- 
;en and sod ium lines are formed primarily in these flows 
Calvet et al. 2000) . A correlation between the emitting 
area of the accretion flow and the magni tude of the mass 
accret ion rate has been suggested by iMuzerolle et al.1 
(2001) as the explanation for the well-established empir- 
ical correlation between infrared hydrogen line luminosi- 
ties and accretion lum i nosities (IMuzerolle et al.l 119981 : 
iFolha fc Emerson! 120011 : iNatta et al.l l2004h . The mod- 
els of hydrogen line formation in magnetospheric flows 
predict that hydrogen line luminosities are primarily de- 
termined by the surface area of the accreting gas, not the 
density in the flow. The suggestion is that objects with 
higher accretion rates require larger emitting areas for 
their magnetospheres than objects with smaller accretion 
rates. Since more extended magnetospheres are expected 
on theoretical grounds in objects with lower disk accre- 
tion rates, a further suggestion is that high-accretion- 
rate objects have wider azimuthal coverage of accreting 
columns. The red absorption profiles of He I A10830 
give new insight into this phenomenon, since we have 
a clear indication of very extended and wide flows in 
stars with low accretion rates. For example, our lim- 
ited phase coverage of the edge-on system AA Tau shows 
that at the same time exceptionally strong red absorp- 
tion at He I A10830 is observed, requiring extensive but 
dilutely filled accretion flows, the hydrogen P7 profile is 
weak, narrow, and symmetric, suggesting a small mag- 
netospheric emitting area if it is formed in the accre- 
tion flow. We anticipate that time-monitoring campaigns 
combining profile monitoring of both He I A10830 and the 
immediately adjacent P7 line of hydrogen will provide a 
definitive assessment of the size and azimuthal coverage 
of the funnel flow and possibly also clarify the origin of 
the correlation between infrared hydrogen line luminosi- 
ties and the accretion luminosity. 

6. CONCLUSIONS 

We have probed the geometry of magnetospheric ac- 
cretion in classical T Tauri stars by modeling red ab- 
sorption at He I A10830 via scattering of the stellar and 



25 



veiling continua. Between 2001 and 2007, we acquired 
81 1-^m spectra of 38 CTTS spanning the full observed 
range of mass accretion rates. Of the 38 stars, 1 of 9 
with ry > 0.5 and 20 of 29 with ry < 0.5 show red 
absorption at He I A10830 that extends below the 1-^m 
continuum in one or more spectra, demonstrating that 
red absorption from magnetospheric accretion is rare in 
objects with high veiling but is found in about two thirds 
of objects with moderate to low veiling. The red absorp- 
tions can be strong, deep, and broad, with equivalent 
widths up to 4.5 A, maximum penetrations into the 1- 
/zm continuum up to 61%, and widths at one quarter of 
the absorption minimum up to 320 km s _1 ; furthermore, 
they tend to be strongest in stars with the lowest veilings. 

We model the red absorption by assuming that an ax- 
isymmetric dipolar accretion flow scatters photons from 
the star and from hot zones in the accretion-heated pho- 
tosphere that produce the l-/xm veiling and have filling 
factor /. Testing a range of magnetosphcre widths and / 
consistent with shock filling factors from the literature, 
we find that about half of the absorption profiles can be 
explained by dipolar flows in which the size of the flow 
is consistent with the size of the shock filling factor /. 
Weak absorptions in stars with weak veiling and inter- 
mediate absorptions in stars with intermediate veiling 
are explained by such flows, but strong absorptions in 
stars with little to no veiling arc not. 

We introduce the concept of dilution as a means of 
producing a strong red absorption while keeping the fill- 
ing factor and thus the veiling low. In a diluted flow, 
the magnetosphere can extend over a wide range of radii, 
with a large covering factor on the stellar surface, but this 
volume is incompletely filled by accreting gas. Instead of 
a single thick flow, we posit multiple nested streamlets 
with a total filling factor small enough for a low veiling, 
but each with an intrinsic thermal or turbulent width 
sufficient to scatter photons as though the entire volume 
were filled, thereby yielding a large red absorption. The 
multiple streamlets can also explain how helium is ion- 
ized through the entire flow, rather than just the skin 
of a thick flow. Large, dilutcly filled accretion flows are 
necessary for about half of the objects, some of which 
require accreting streamlets to connect to the disk over a 
range from 2 R* out to or beyond corotation. A few stars 
show such deep absorption at redward velocities exceed- 
ing 50% of the stellar escape velocity that flows near the 
star with less curvature than a dipolar trajectory seem 
to be required. 

The frequency of He I A10830 red absorption is also in- 
formative. Our limited temporal coverage suggests that 



the frequency of helium absorption differs in stars with 
high and low veiling. Red absorption at He I A10830 is 
far more common in stars with low veiling. When it is 
absent from these stars, it is sometimes because helium 
emission from another source such as a wind fills it in 
and sometimes because of inhomogeneous azimuthal cov- 
erage of accreting magnetic columns. Among stars with 
high veiling (ry > 0.5), red absorption at He I A10830 
is rarely seen. If these stars had accretion geometries 
similar to those of the low-veiling stars, they would be 
expected to have extremely strong red absorptions. Even 
if the absorption were filled in by emission from the ac- 
cretion flow, the stars would still be expected to show red 
absorption at high velocities. In the high-veiling stars, 
the paucity of He I A10830 red absorption, the presence 
of He I A 10830 emission and blue absorption that sug- 
gest formation in accretion-powered stellar winds, and 
the weakness or absence of narrow-component He I A5876 
emission from an accretion shock lead us to suggest that 
the magnetospheric accretion structure may be crunched 
or otherwise reduced in CTTS with the highest disk ac- 
cretion rates. 

We find the study of He I A10830 red absorption due 
to infalling gas projected in front of the star to be com- 
plementary to studies of emission lines modeled as aris- 
ing over the full size of the accretion flow. The prox- 
imity of He I A10830 and P7 offer an excellent pair of 
lines for deeper investigation of magnetospheric geome- 
tries through intensive time-monitoring programs that 
can track non-aziumuthal structures as stars rotate. Our 
limited phase coverage of AA Tau demonstrates that this 
approach will be very effective, particularly when cou- 
pled with radiative transfer models that can constrain 
formation conditions for both lines simultaneously. 
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